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SECTION  1 


Introduction 


In  recent  years,  electromagnetic  scattering  from  complex  structures  has  become  an 
important  field  of  research.  Accurate  and  efficient  solutions  to  scattering  from  com¬ 
plex  structures  can  enhance  the  ability  to  understand  how  targets  can  be  detected 
in  the  presence  of  clutter.  Many  approaches  are  available,  but  for  high  frequency 
scattering,  the  Uniform  Theory  of  Diffraction  (UTD)  is  one  of  the  most  efficient 
and  accurate  techniques.  UTD  can  provide  the  characteristic  scattering  properties 
of  local  features  of  the  target  such  as  frequency  response,  polarization,  and  angular 
behavior. 

It  is  well  known  that  due  to  the  highly  local  nature  of  high  frequency  scattering, 
the  solution  to  complex  structures  can  be  modeled  approximately  as  a  superposition 
of  scattered  fields  from  simpler  shapes  such  as  plates,  ellipsoids,  cylinders,  dihedrals, 
trihedrals,  et  cetera.  These  can  be  used  to  represent  component  parts  of  vehicles 
such  as  trucks,  aircraft,  etc.,  as  well  as  natural  objects  such  as  tree  trunks  and 
branches. 

The  object  of  this  research  is  to  find  far  zone  backscatter  and  bistatic  scattered 
fields  from  some  of  the  most  common  shapes  which  require  second  and  higher  order 
interactions  to  get  an  accurate  result  over  a  wide  range  of  angles  and  frequencies. 
The  basic  shapes  important  to  this  discussion  are  the  dihedral,  bruderhedral,  tri¬ 
hedral,  top-hat,  bi-cone,  parallel  plate  cavity,  cake  pan  (shallow  circular  cavity). 
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doniit  and  miilti-joined  cylinders.  All  the  shapes  considered  are  assumed  to  he  per¬ 
fectly  conducting  except  the  multi-joined  cylinders  which  may  be  made  of  dielectric 
materials.  For  these  shapes,  the  most  important  higher  order  interactions  will  be 
determined.  Where  possible  the  expressions  will  be  simplified  to  show  the  connec¬ 
tion  between  the  UTD  results  and  the  more  classical  PO  approximations.  It  should 
be  pointed  out  that  many  of  the  targets  do  not  have  easily  found  references  to  their 
scattering  properties. 

Some  background  information  on  UTD  (1]  is  presented  in  Section  2.  The  double 
diffraction  formulation  presented  here  is  derived  in  Reference  [2].  In  Section  3,  the 
UTD  analysis  of  a  two-dimensional  dihedral  is  discussed.  ParaDel  ray  formulation  of 
the  UTD  is  used  here  for  the  far  zone  problem.  This  means,  in  two  dimensions,  that 
only  edge  diffractions,  their  images,  and  various  combinations  of  double  diffractions 
are  needed.  This  is  a  very  efRcient  and  accurate  means  of  analyzing  the  far  zone 
problem.  It  avoids  any  integrations  and  non  physical  “imposed”  edged  conditions. 
Basic  dihedral  mechanisms  are  studied,  and  it  is  shown  how  the  various  interac¬ 
tions  ensure  continuity  of  the  total  field.  It  is  found  that  excellent  comparisons  with 
Method  of  Moments  results  are  obtained  using  fields  of  up  to  third  order  interactions 
for  dihedrals  of  60°  or  greater.  Also,  it  is  observed  that  this  solution  is  good  for  di¬ 
hedrals  that  are  as  small  as  about  one  wavelength  (maximum  width  of  the  dihedral). 
Higher  order  interactions  are  required  to  get  results  for  less  than  60°  dihedrals  or  for 
smaller  sizes.  For  the  case  of  maximum  scattering  from  a  general  angled  dihedral, 
the  UTD  terms  that  are  dominant  in  between  the  strips  of  the  dihedral  are  reduced 
to  simple  forms. 

UTD,  being  a  ray  technique,  predicts  infinite  fields  at  caustics.  It  is  therefore 
iiecessary  to  introduce  caustic  corrections  to  obtain  a  bounded  field  at  the  caustics. 
Caustic  corrections  of  first,  second  and  higher  order  are  discussed  in  Section  4.  In 
Section  5  and  6,  the  two-dimensional  dihedral  solution,  modified  by  the  spread  factor 
and  the  caustic  corrections  discussed  in  Section  4,  is  used  to  get  results  for  the  top- 
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liat  and  Id-rone,  respectively.  It  is  shown  that  the  nirvatnre  modifications  and  the 
caustic  corrections  to  the  two-dimensional  solution  gives  appropriate  results  when 
compared  with  the  body  of  revolution  results  even  for  very  small  sizes. 

A  modification  to  the  two-dimensional  dihedral  is  the  parallel  plate  cavity.  The 
parallel  plate  cavity  modified  for  curvature  and  caustics  effects  is  the  cake  pan  (shal¬ 
low  circular  cavity).  In  Section  7,  the  parallel  plate  cavity  and  the  cake  pan  are 
discussed.  For  the  cake  pan,  ray  terms  which  contribute  only  in  the  principal  plane 
are  included.  It  is  found  that  for  shallow  and  wide  circular  cavities  good  results 
can  be  obtained,  over  a  wide  range  of  angles,  by  including  only  the  ray  terms  that 
contribute  in  the  principal  plane. 

For  the  donut,  in  Section  8,  the  first  and  second  order  GO  terms  and  its  caustic 
corrections  are  discussed.  No  creeping  wave  terms  are  included  in  this  solution. 

In  Section  9,  techniques  developed  for  the  two-dimensional  dihedral  solution  and 
the  cylinder  solution  is  also  used  to  get  backscatter  results  from  multi-joined  cylin¬ 
ders.  This  model  is  used  to  investigate  the  scattering  from  geometries  similar  to  tree 
trunks  and  branches.  In  particular,  two  types  of  “tree”  crook  models  are  examined. 
The  UTD  results  are  compared  with  the  measurements  and  reasonable  agreement  is 
found. 
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SECTION  2 


Theoretical  Background 

1  Introduction 

This  chapter  contains  a  brief  description  of  the  techniques  of  Geometrical  Optics 
(GO),  and  the  Uniform  Theory  of  Diffraction  (UTD)  needed  for  the  solution  of  the 
configurations  analyzed  in  this  work.  Details  on  these  high  frequency  approximate- 
solution  methods  can  be  found  in  References  [1,  3,  4,  5].  Also,  a  recently  developed 
far- zone  double  diffraction  coefficient  is  presented  [2]. 

It  is  assumed  that  the  fields  exhibit  a  time  dependence  of  e-'” '  which  is  suppressed 
throughout  this  work. 

2  Geometrical  Optics  (GO) 

A  brief  summary  of  geometrical  optics  is  given  here.  The  propagation  of  electromag¬ 
netic  energy  at  high  frequencies  through  isotropic,  lossless  media  can  be  described 
using  GO.  Thus  GO  provides  a  high  frequency  approximate  solution  to  the  incident, 
reflected  and  refracted  fields.  In  our  problem,  we  will  be  concerned  with  the  reflected 
field.  The  GO  can  be  obtained  from  an  asymptotic  solution  of  Maxwell’s  equations 
and  it  corresponds  to  the  leading  term  of  the  Luneberg- Kline  asymptotic  expansion 
for  large  values  of  angular  frequency  [1).  According  to  GO,  electromagnetic  energy 
can  be  viewed  as  travelling  along  well-defined  paths  known  as  rays.  The  ray  path,  in 
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any  continuous  mediiim,  may  be  determined  using  Fermat’s  principle  wbicli  states 
the  energy  will  flow  along  the  path  of  shortest  electrical  length  between  any  two 
points.  The  shortest  electrical  path  is  one  which  results  in  the  shortest  propagation 
time  between  two  points.  From  Fermat’s  principle,  the  law  of  reflection  and  the  law 
of  refraction  can  be  obtained.  In  geometrical  optics,  phase  is  proportional  to  the 
distance  along  the  ray  path  from  some  reference  point,  and  the  amplitude  is  gov¬ 
erned  by  the  conservation  of  power  in  a  tube  of  rays.  Figure  1  shows  one  such  tube 
of  rays.  The  distances  p\  and  pa  are  from  a  fixed  reference  point  to  the  “caustics”, 
where  a  family  of  rays  merge  or  intersect. 

In  geometrical  optics,  like  all  other  ray  techniques,  the  fields  can  not  be  evaluated 
at  the  caustics.  At  caustics,  the  approximation  of  power  conservation  in  a  tube  of 
rays  is  no  longer  valid.  The  field  near  a  caustic  appears  to  become  infinite  as  a  finite 
amount  of  power  gets  squeezed  into  a  vanishing  area.  The  fields  near  caustics  must 
be  found  by  other  means.  Chapter  4  discusses  some  of  the  caustics  encountered  in 
this  work. 


2.1  GO  Reflected  Field 


The  GO  reflected  field  from  a  smooth  conducting  convex  surface  with  radius  of 
curvature  greater  than  A,  see  Figure  2,  is  given  as 


t{f)  =  E'{Qr) 


R\ h— 


P]P2 


(1) 


where 

Qr  =  reflection  point, 

(Qr)  =  incident  field  at  Qr, 
s’’=  distance  from  Qr  to  receiver, 

R  =  the  dyadic  reflection  coefficient, 
R  =  R„e\e\  -f  jRhef|e[j, 

R.,h  =  Tl, 
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REFERENCE  POINT  FIELD  POINT 

RAY 

Figure  1:  An  astigmatic  tube  of  rays. 
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fi*UNIT  OUTWARD  NORMAL  VECTOR  TO 
THE  CONVEX  SURFACE  AT  0. 

II  JL 

l[l  « 


Figure  2:  Reflection  off  a  convex  conducting  surface. 
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Pi .2—  principal  raHii  of  curvature  of  reflected  wavefront,  and 


J_  =  1  fl  .  ^  [(022^  +  (0.2^  (02.f  +  (0..f' 

Pl.2  2  {p\  ^  pU  ^  (0|2  [  R,  R2  . 

.  1  I  /I  _  ±y  4cOsg»  [(022^  +  (0|2f  (02.f  +  (01lf' 

2\vp‘.  pi/  vpi  Pi/  I0i"  I  i2,  i22  , 

.4C0S*d'  r/(022f  +  (0.2f  .  (02.f  +  (0.1f^'  4101*11  ^ 

10|-  [V  R2  /  ~  RiR2\]  ’ 

where  the  plus  sign  is  associated  with  p\  and  the  minus  sign  with  pji 

0^*  =  X).U,  (3) 


or 


i',.172 
X2.U^  X2.V2 


(4) 


JVj  and  A'2  are  the  principal  directions  of  the  incident  wavefront  on  the  curved 
surface  S  at  Qr  with  the  principal  radii  of  curvature  p\  and  p\^  and  O'  is  the  angle 
of  incidence.  (See  Figure  3)  U\  and  U2  are  unit  vectors  in  the  principal  directions  of 
S  at  Qr  with  the  principed  radii  of  curvature  R\  and  R2. 

For  the  far-zone  case,  the  Equation  (2)  can  be  simplified  as 


1  _  1  sin*  O2  sin*  9\ 

p\  2  COSP’  Rx^R2 


where  9\  is  the  angle  between  the  direction  of  the  incident  ray  i'  and  {/],  and  62  is 
the  angle  between  3*  and  U2- 


3  Uniform  Theory  of  Diffraction  (UTD) 

Keller’s  generalization  of  fermat’s  principle  allows  one  to  include  a  class  of  rays 
diffracted  from  electrical  and/or  geometrical  discontinuities  on  a  scattering  or  a 
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Figure  3:  Geometry  of  the  description  of  the  wavefront  reflected  from  the  curved 
surface  S. 
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radiating  object;  these  diffracted  rays  exist  in  addition  to  the  rays  of  GO  which 
describe  incident,  reflected  and  refracted  rays. 

In  GTD,  the  initial  value  of  the  diffracted  rays  are  obtained  by  multiplying  the 
incident  field  with  a  diffraction  coefficient  at  the  point  of  diffraction,  analogous  to  the 
reflection  coefficient  for  the  reflected  ray.  The  expression  for  the  diffraction  coefficient 
depends  on  the  scattering  mechanism  (edge  diffraction,  tip  diffraction,  etc.).  The 
application  of  GTD  is  limited  only  by  the  availability  of  accurate  coeflUcients.  The 
Uniform  Theory  of  Diffraction  (UTD)  is  an  improved  version  of  the  GTD  in  that 
the  diffraction  coefficients  remain  valid  at  and  near  the  shadow  boundaries,  where 
the  Keller’s  original  coefficient  had  failed. 

Two  major  contributors  to  the  diffracted  field  considered  by  UTD  are  diffrac¬ 
tion  from  the  edge  and  the  diffraction  from  a  curved  surface  (or  the  creeping  wave 
diffraction).  In  this  work,  only  the  edge  diffraction  will  be  considered. 


4  UTD  Edge  Diffracted  Field 

The  singly  edge  diffracted  UTD  field  as  developed  by  Kouyoumjian  and  Pathak  [1] 
is  now  introduced.  The  three-dimensional  wedge  diffraction  geometry  is  shown  in 
Figure  4.  Let  an  arbitrarily  polarized  field  E'{8')  which  may  be  plane,  cylindrical  or 
spherical  wave  be  incident  on  the  wedge.  Such  a  field  £  {9')  can  be  represented  as 
a  GO  incident  ray  field.  When  an  incident  ray  strikes  the  edge  at  Qr  (say),  then  it 
generates  a  cone  of  edge  diffracted  rays  which  emanate  from  so  that  the  cone  half 
angle  =  angle  made  be  the  incident  ray  into  the  edge  tangent  at  Qg  (i.e.,  /?„  = 

This  is  a  consequence  of  the  generalization  of  the  fermat’s  principle  to  the  incident 
edge  diffracted. 

The  authors  of  [1]  expressed  the  dyadic  diffraction  coefficients  in  a  ray  fixed 
coordinate  system  which  is  centered  at  the  diffraction  point  Qg.  The  orthogonal 
unit  vectors  associated  with  the  ray  fixed  coordinate  system  are  defined  as: 


•  / 
s 


(6) 
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Figure  4:  Geometry  for  the  three-dimensional  wedge  diffraction  problem. 


and 

P„y.  (7) 

These  vectors  are  shown  in  Figure  4.  The  diffracted  field  from  a  semi-infinite 
perfectly-conducting  curved  wedge  is  given  by  the  following  expression: 

where 

Qf.  =  point  of  diffraction  on  the  edge, 

E' (Qr)=  incident  field  at 

s  =  distance  from  Qr  to  receiver, 

p**  =  caustic  distance  for  the  diffracted  ray, 

^  =  the  dyadic  edge  diffraction  coefficient, 

(9) 

or  in  matrix  form 


The  ray  fixed  coordinate  system  ensures  that  the  diffraction  coefficient  dyadic,  D, 
matrix  is  diagonal. 

The  soft  diffraction  coefficient,  D.,,  is  associated  with  the  component  of  E-field 
that  is  parallel  to  the  edge,  i.e., 

zE  =  0.  (11) 

The  hard  diffraction  coefficient,  £?/,,  is  associated  with  the  component  of  E-field  that 
is  perpendicular  to  the  edge,  i.e., 

hv(^-ff)  =  0.  (12) 
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The  soft  and  hard  diffraction  coeffirients  are  defined  as 


Dt,h  — 


where 


2n\/27rfc  sin  fio 


+  cot  2n~^  ^ 


T 

4-  cot 


I.  ~  (2n-irN^  —  (f3)\ 

a±(/9)  =  2cos=  I - J— 

and  are  integers  which  most  nearly  satisfy  the  equations 


(13) 


(14) 


2-n:nN'*'  —  (/3)  =  ir, 


(15) 


and 

2irniV-  -  {0)  =  -w,  (16) 

with 

0  =  0^  =  (17) 

F(x),  where  x=kLa,  is  a  transition  function  defined  as 

/’(x)  =  2j\y/i\  e"-'^  dr  (18) 

The  magnitude  and  phase  of  F(x)  is  shown  in  the  Figure  5.  For  x  small  (X  <  0.3), 

F(i)  ~  —  2xe^*  —  (19) 

When  0.3  <  x  <  3.5,  linear  interpolation  is  used  as  an  efficient  numerical  means  of 
determining  F(x)  so  that 


F(x)  =  F(x„)  +  4„(x  -  In). 


(20) 
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Figure  5:  Magnitude  and  phase  of  F(x). 


For  X  >  5.5  one  finds  that 


r./  X  ,  J  3  .15  75 

F(x)  ~  1  H - - —  7 - h 

^  v-;  ^  o_ 


2x  4z''*  8x  ‘  16x ' 


The  diffracted  caustic  distance,  p’^,  is  defined  as 

1  In,  -  (s'  -  s) 

-  \  =■  \  (22 

p^  p\.  a,  sin^ 

L  is  a  distance  parameter,  where 

t"'  = »  ^ 

pt  (P,  +5j(P2  +«j 

The  parameters  p'/  are  the  (incident,  reflected)  radii  of  curvature  at  Q,  in  the  plane 
containing  the  (incident,  reflected)  ray  and  e  is  the  unit  vector  tangent  to  the  edge  at 
Qff  Tif  is  the  associated  unit  normal  vector  to  the  edge  directed  away  from  the  center 
of  curvature,  a,  >  0  is  the  radius  of  curvature  of  the  edge  at  Qr  and  is  the  angle 
between  the  incident  ray  and  the  tangent  to  the  edge  as  shown  in  Figure  4.  The 
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unit,  vectors  f'  and  a  are  in  the  direction  of  incidence  and  diffraction,  respectively. 
p’,’2  are  the  principal  radii  of  curvature  of  the  (incident,  reflected)  wavefront  at  Q,.. 
The  superscript  n,  o  on  L  denote  that  the  radii  of  curvature  are  determined  at  the 
reflection  shadow  boundaries  of  the  n  and  o  faces,  respectively.  The  reflected  field 
caustic  distance  in  the  plane  of  the  reflected  ray  and  edge  tangent  is  given  by 


Pc 


1  2(n.nc)  (i'.n) 

P'c 


(24) 


ar  sin*  /?„  ’ 

where 

n  =  normal  to  the  surface  at 

hr  —  normal  to  the  edge  curvature  at  Qr, 

s’  —  incident  ray  direction, 

a,.  =  radius  of  the  edge  curvature. 

The  diffraction  coefficient  in  Equation  (13)  may  be  simplified  for  the  special  case 
when  the  source  aad  the  receiver  are  in  the  far  field  so  that  F(x]  is  approximately 
equal  to  1.  For  this  case  the  diffraction  coefficients,  Z),./,  reduce  to  the  Keller’s 
diffraction  coefficients  in  Reference  [3],  or 

[cot  (^)  +  cot 


{„t(^),„t(l^)}].  ,.5) 


which  may  also  be  written  as 


— e  sin  - 

_ _ 

n-\/2irfc  sin  Po 


1 

“T 

1 

;)  COs(= 

(26) 


The  diffraction  coefficient  may  also  be  simplified  for  the  special  case  when  n=2 
(half-plane) 


—  e 


D,,h{L,<f>,<f>\P„)  =  y-=Y 


2\2Trk  sin 


F\kLa{0~)]  ^  F[kLa{p^)] 


cos(^-/2)  cos(i9+/2) 
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(27) 


5  Scattering  from  a  Strip  in  the  Far-zone 


A  study  of  the  scattering  from  a  strip  is  considered  because  the  strip  is  the  basic 
building  block  for  the  structures  discussed  in  this  work.  For  the  case  of  far-zone  illu¬ 
mination  and  observer,  the  first  order  bistatic  scattered  field  from  the  strip  contains 
only  the  diffracted  field  from  two  edges.  The  GO  information  need  not  be  specifi¬ 
cally  included,  since  it  is  contained  in  the  GTD  far  zone  diffracted  fields.  Consider 
the  geometry  shown  in  Figure  6.  The  length  of  the  strip  is  given  as  1.  The  basic 
geometrical  parameters  needed  are 


<f>i  -  <f>\  =  <f>  -  (f>' 


4>\  +  4>\  =  <l>  + 

rr  +  <t>„-  <f>'  ,  0  <  (p'  <  T 

2ir  <l>„  -  <f>'  ,  IT  +  4>.>  <  (p'  <2ir 


\  IT  +  <p„-  <p  ,0  <  <p  <  n  (po 
<p2  =  < 

y  Zn  <p„  -  <p  ,n  +  <p„<<p<2n 

Ki.  =  ‘V’,!, 


(28) 

(29) 

(30) 

(31) 

(32) 


where 


and 


I, 


E„  ,  E  parallel  to  wedge 

< 

H„  ,  E  perpendicular  to  wedge 


(33) 


ui  =  u\  , _ 

^  ^  2v27rfc 


+ 


cos 


U\- 


1 


(^) 


i-T^j  cos 

+« "!*  ) 


.  cos 


COS 


(34) 


* 

It  is  observed  that  both  diffraction  coefficients  become  infinite  at  the  reflection 
shadow  boundaries  (i.e.  <p\  <p\  —  w  and  <^2  +  (^2  =  tt;  or  -f  ^'  =  2<p„  +  tt)  and  at 

the  incident  shadow  boundary  (forward  scatter)  (i.e.  <p\  ~  (p\  —  ir  and  <p2  —  (p'^  =  tt). 
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Figure  6:  Geometry  of  the  strip. 
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Using  Equations  (28)  through  (31)  into  Equation  (34)  and  constraining  to  the 
region  0  <  ^  <  ir  +  ^„  and  0  <  <^'  <  ir  -f  gives 

...  ... -e-J*/'  /  1  1  \ 


Ui  =  U\ - = 

h  h  2VS 


Trk  Vcos(^] 


2vCT  \cos  (^)  ^  cos  j  • 

Notice  that  for  this  range  the  scattering  along  the  incident  shadow  boundary  is  not 
possible.  Combining  terms  results 

h  h  2\/2iilfc 


2  j  sin  1 

j- 

kl cos  ^ 

cos  1 

(^)) 

COS 

) 

2 cos  (kl cos  cos  ) 

In  general,  the  case  of  infinity  along  the  reflection  shadow  is  more  common  and 
is  treated  in  detail  to  show  how  a  continuous  result  is  obtained.  Near  the  reflection 
shadow  boundary 


-f  —  2<^„'j  _ 


cos(Tr/2  ±  e)  ~  q^c. 


Using  Equation  (37)  in  Equation  (36)  gives 

l/i  =  (ijfclcos  + - 7X^1  • 

f  /.  \  ^  \  2  )  cos(^)j 


For  broadside  backscatter  case 

Ui  =  ui.  —7=  [±jki  + 1]  =  t/r  +  u'-  (39) 

where  UH”  the  PO  term  and  U'^  is  the  correction  to  PO.  Simplifying  the  U^'  term 

h  <  h 


ur  =  ^u\--=ki 

h  v27rfc 
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_ 


2 


(41) 


gives  the  expected  PO  result  for  backscatter  from  broadside  of  a  2-D  strip. 

Results  of  fields  backscattered  from  a  strip  are  presented  in  Figures  7  and  8.  The 
UTD  solution  is  obtained  by  the  above  formulation.  Comparing  the  UTD  result 
with  the  moment  method  result,  it  is  observed  that  using  the  single  edge  diffracted 
fields  from  the  edges  gives  a  good  approximation  to  the  exact  solution  for  the  main 
and  a  few  side  lobes.  For  the  T case,  we  note  that  near  <^  =  0"  and  180"  the  first 
order  UTD  solution  is  not  accurate.  This  is  a  consequence  of  the  higher  order  terms 
not  being  taken  into  account.  However  for  the  TMz  case  the  higher  order  terms  are 
negligible.  See  References  [6,  7,  8,  9]  for  higher  order  effects  on  the  strip. 

6  Double  Diffraction  Formulation 

Previous  solutions  of  the  field  diffracted  by  double-edge  structures  used  the  UTD 
in  its  original  formulation  argumented  by  slope  diffraction  (especially  for  the  soft 
polarization)  to  give  useful  results  when  the  distance  between  the  two  edges  is  very 
large.  However,  difficulties  are  present  whenever  the  second  edge  is  located  in  the 
transition  region  of  the  first  edge  and  the  diffracted  field  is  calculated  in  the  transition 
regions  adjacent  to  the  incident  or  reflection  shadow  boundaries  of  the  second  edge. 
These  problems  are  caused  by  the  rapid,  spatial  variation  together  with  a  non-ray 
optical  behavior  of  the  field  incident  on  the  second  edge  after  diffraction  from  the 
first.  It  is  obvious  that  the  angular  range  of  the  overlapping  transition  regions 
broadens  as  the  distance  betwen  '.e  edges  decreases.  Part  of  the  difficulties  are 
removed  by  a  spectral  extension  of  the  UTD  in  References  [6]  and  [7).  This  solution 
is  uniformly  valid  at  any  incidence  aspect,  but  it  is  restricted  to  those  observation 
aspects  which  coincide  with  diffraction  shadow  boundary  (DSB)  or  the  diffraction 
reflection  shadow  boundary  (DRSB).  (See  Figure  25.) 
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A  closed  form  solution  for  the  doubly  diffracted  field  in  the  far  zone  from  a  pair  of 
parallel  wedges  is  presented  in  Reference  [2].  As  shown  in  Section  3,  these  expressions 
exhibit  the  proper  discontinuities  and  singularities  so  that  when  combined  with  the 
singly  diffracted  field,  the  field  is  continuous  and  uniform.  For  joined  wedges,  shown 
in  Figure  9,  the  angle  parameters  are  defined  as 


f  -  ^’i); 


(42) 


and 

=  IT  +  (-l)''(n2ir  -  (^2);  g  =  1,2.  (43) 

The  pattern  function  of  the  field  doubly  diffracted  from  Qj  and  Q2  is  given  as 


2  2 


(44) 


/)=l  g=l 


By  only  retaining  the  first  nonvanishing  term  in  the  asymptotic  expansion  for  kd 
large,  the  following  closed  form  expressions  are  obtained 

(-i)''(-i)’ 


P.(*r.^) 


CLpdq 


F[kda,.]  -  F{kda,} 


Op  o. 


(45) 


— 

)  cot  { - I 

\2nj  j 

>  \2n2J 

(46) 


for  a  TM  plane  wave  (with  the  electric  field  parallel).  For  a  TE  plane  wave  (  with 
the  electric  field  perpendicular), 

F,.($p,$,)~-— i - 

47ryn,n2 

agF[kdaj,]  -  apF[fcda,]  / 

flp  \  ^ kd  / 

and 

o,  =  2cos*|($,  —  TT  —  (  — l)'2Ti7rlV,)/2]  (47) 

where  n  =  ni  for  i  =  p,  n  =  712  for  i  =  q,  and  N,  is  the  integer  that  most  nearly 
satisfies 

V 


JV,  =  (-!)' 


2mr 


(48) 
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F(x)  is  the  transition  function  that  has  been  defined  for  the  single  diffraction  case. 

If  the  two  wedges  are  separate,  the  above  equations  are  slightly  varied.  For  the 
separated  wedges  shown  in  Figure  10,  the  angle  parameters  are  defined  as 

IT  +  p=1,4,  (49) 

=  IT  +  (-1)'»(^2  -  9=1>4,  (50) 

with  q,  =  1  for  1=1,  2  and  ijj  =  —  1  for  i=3,  4.  The  pattern  function  of  the  field 
double  diffracted  from  Qi  and  Q2  is  given  as 

Pi  ~  E  E  ^  |P.  (»,..♦,)  +  Pi.  (^r.*,)] ,  (51) 

;.=  U=I  ^ 

where  x«  =  Vi  TM  case  and  Xt  =  1  in  the  TE  case. 

7  Limit  Forms  of  Double  Diffraction  Coefficients 

When  a,,  and  a,  are  such  that  they  are  approximately  equal  but  not  close  to  zero, 
Equations  (45)  and  (46)  tend  to  become  of  indeterminate  form.  This  can  cause 
numerical  problems  on  the  computer.  Equations  (52)  and  (53)  are  parts  of  Equa¬ 
tions  (45)  and  (46),  respectively,  that  may  have  such  a  behavior.  To  get  around  this 
problem  on  the  computer,  limit  forms  are  used  in  the  approximate  region.  Defining 


y  _  F{kda,]  -  F[kdaq] 

(62) 

dp  dq 

and 

a^F\kdaj]  -  apFffcda,] 

—  » 

(53) 

dp  —  dq 

and  noting  that 

Wlim  + 'i  -  =  kiF'\Ma]  U.. 

kdaa  +  t-kdaa  1 

(54) 

where  the  /  denotes  diffraction  with  respect  to  the  argument  of  the  transition  func¬ 
tion. 

Ip  Ip  "  -/(P)>  (^5) 
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Figure  9:  Joined  wedges  geometry. 
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and 


F[x]  =  2jy/ie^^  /  dr,  (56) 

J  y/i 

equations  for  T,  and  T/,,  that  are  valid  when  a,,  approaches  a,  may  be  derived.  From 
Equations  (54)  to  (56)  one  obtains, 

dF 


dx  ^ 


1  .  9 

1  -  (-4=  +  2jy/i)e^^  /  dr 

\  X  J  y/x 


^  [2jx  -  (1  +  2jx)  F{x)] , 


and  as  Op  — >  a,,  from  Equation  (52) 

lim  T,  =  1 1  — I-  jAjd )  F\kdaj]  —  jkd 

"<»— "p  L\2“p  / 


Equation  (53)  may  be  rearranged  as 


_  /'[A-rfool 


Th  —  Ofndq 


dp  dq 


and 

lim  Th  =  {kd^apdp— 
a,-<ip  dx 

Using  Equation  (58)  in  Equation  (61) 


lim  Th  =  kdd, 

a,— Op 


(•^  2kdd^ 


(57) 

(58) 

(59) 


“I 

> 

(60) 

[fN] 

(61) 

X 

^  F(i)  -  j 

■ 

(62) 

Thus,  when  dp  and  a,  are  almost  equal  but  not  close  to  zero.  Equations  (59)  and 
(62)  may  be  replaced  by  Equations  (45)  and  (46)  for  T,  and  Th  to  insure  correct 
results  and  prevent  computer  errors. 
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SECTION  3 


UTD  2-D  Analysis  of  a  Dihedral 

1  Introduction 

The  dihedral  has  been  in  use  as  a  scattering  device  for  many  years.  Some  studies 
conducted  on  the  dihedral  at  The  Ohio  State  University  ElectroScience  Laboratory 
are  given  in  References  [10,  11,  12,  13].  They  entail  mostly  measurement  studies 
with  some  basic  analyses.  More  recently,  renewed  interest  in  the  dihedral  has  led 
to  more  detailed  analytic  studies.  A  Physical  Optics  for  a  right  angled  dihedral  is 
analyzed  with  full  polarization  scattering  matrix  in  Reference  [14] .  Physical  Optics 
and  the  Physical  Theory  of  Diffraction  are  used  to  predict  the  results  for  dihedrals 
in  Reference  [15].  In  Reference  [16],  a  near  zone  type  UTD  plus  an  “imposed”  edge 
diffraction  extension  is  used  to  get  the  backscatter  results  for  the  dihedral  in  the  “far 
zone”  by  extending  the  distance  parameter  for  the  near  zone  case.  The  use  of  UTD 
to  study  plate  interactions  in  the  far  zone  is  conducted  in  Reference  [17].  This  type 
of  solution  is  improved  and  extended  to  specific  dihedral  problems  in  Reference  [13], 
and  also  briefly  discussed  in  this  report  as  a  foundation  for  the  other  scatter  objects 
analyzed  in  this  work. 

As  seen  in  the  previous  section,  the  field  of  a  strip  can  be  obtained  by  including 
the  diffracted  fields  from  the  two  edges,  and  the  infinity  from  edge  diffractions  com¬ 
bines  to  give  a  bounded  result  in  the  specular  direction.  When  more  than  one  strip 
is  used  in  a  model,  the  singly  diffracted  field  from  the  edges  will  act  as  a  near-field 
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source  which  can  create  reflected  and  diffracted  fields  from  the  other  strip.  Thus, 
interactions  between  the  strips  are  necessary  to  provide  a  more  accurate  pattern.  In 
general,  these  interactions  are  present  only  over  a  limited  region  which  may  be  deter¬ 
mined  geometrically.  These  higher  order  interactions  also  provide  the  cancellation 
of  infinities  to  provide  a  bounded  result.  It  will  be  shown  that  interactions  up  to 
third  order  are  sufficient  to  obtain  results  with  engineering  accuracy  for  the  case  of 
two-dimensional  dihedrals  of  60°  or  greater  and  for  sizes  larger  than  one  wavelength 
in  total  length. 

The  UTD  interactions  that  have  been  included  in  the  analysis  which  are  contained 
in  the  2-D  computer  code  are  shown  in  Figure  11.  The  single  interaction  fields  in¬ 
cluded  are  simply  edge  diffraction  (D)  terms.  The  second  order  interaction  fields  are 
the  reflected  -  diffracted  (RD),  diffracted  -  reflected  (DR),  and  diffracted  -  diffracted 
(DD)  (double  diffracted)  terms.  The  third  order  interaction  fields  are  the  reflected 
-  reflected  -  diffracted  (RRD),  diffracted  -  reflected  -  reflected  (DRR),  reflected  - 
diffracted  -  reflected  (RDR),  reflected  -diffracted  -  diffracted  (RDD),  diffracted  - 
diffracted  -  reflected  (DDR),  and  diffracted  -  reflected  -  diffracted  (DRD)  terms. 
Triple  diffractions  have  not  been  included.  Any  of  these  terms  could  be  critical  es¬ 
pecially  at  a  shadow  boundary.  As  a  general  rule,  by  knowing  the  position  of  the 
source  and  the  receiver,  and  by  drawing  ray  diagrams,  it  is  simple  to  predict  the 
important  terms  by  observing  the  directions  of  the  bounces. 

The  solution  discussed  here  is  based  on  a  study  of  strip  interactions  given  in 
Reference  [17]  and  a  study  of  double  diffraction  in  Reference  [2].  The  strip  interaction 
in  the  previous  study,  however,  did  not  provide  the  necessary  accuracy  for  dihedral 
configurations.  The  two  dimensional  geometry  studied  is  shown  in  Figure  12.  The 
characteristic  features  of  our  analysis  allows  that  :  (i)  the  length  of  each  strip  may 
be  different,  (ii)  the  angle  between  the  strips  may  be  varied,  (iii)  back  and  bistatic 
scattering  may  be  treated,  (iv)  linear  and  elliptical  polarization  can  be  handled,  (v) 
the  solution  is  valid  for  small  sizes  in  terms  of  a  wavelength. 


Figure  11:  The  interactions  considered. 
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2  Basic  Dihedral  UTD  Mechanisms 


This  section  discusses  the  basic  UTD  terms  that  provide  the  characteristic  dihedral 
pattern.  In  the  case  of  the  main  beam  of  the  pattern,  the  RD  and  the  DR  UTD 
terms  along  with  the  single  diffracted  (D)  term  from  the  wedge  formed  at  the  junction 
between  the  two  intersecting  strips  provide  the  necessary  fields.  The  specular  from 
the  individual  strip  faces  is  produced  by  single  diffraction  (D)  from  the  ends  of  the 
strips  and  the  DR,  RD  and  DRD  terms  between  the  strips.  The  double  diffraction 
(DD)  across  the  face  of  the  dihedral  can  also  be  important. 

As  mentioned  above,  the  RD  and  DR  terms  along  with  the  D  term  from  the 
junction  provide  the  starting  point  for  determining  the  scattering  from  an  arbitrar¬ 
ily  angled  dihedral.  These  terms,  in  general,  look  quite  complicated.  In  some  special 
cases,  they  can  be  reduced  to  simpler  formulas.  As  expected,  in  the  case  of  backscat- 
ter  for  a  wedge  angle  of  a  =  90°,  the  well  known  PO  result  for  a  dihedral  can  be 
derived. 

From  the  equations  in  Section  2  and  for  angles  derived  from  Figure  13  for  bistatic 
scattering  from  a  general  angled  dihedral,  the  diffracted  term  from  the  junction  is 


u'^  =  u: 


_e-Wi  ( 

<  cot 


*’*  2n\/2irk  \ 
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^TT  -f 
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-f  cot 
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In  Figure  13,  D\R  and  RDi  terms  are  due  to  reflection  from  strip  2.  The  terms 
D2R  and  RD2  are  due  to  reflection  from  strip  1.  The  basic  geometrical  parameters 
for  the  Di  R  term  are 


—  IT  —  (f)  ]  —  TT  —  2Q[  —  ijf 

(64) 

=  27r  —  2q  <f>  —  (f> 

(65) 

^1  ~  ~  — 2ci  (f>  -{■  <f> . 

(66) 
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For  the  RD\  term,  the  angles  are  the  same  as  for  D\R  with  </>  and  <f>'  interrhanged. 
The  term  RD2  has  the  following  parameters 


^2  =  IT  -  (a  +  4>')  ;  =  TT  -  o  + 


(67) 


02  +  ^2  ~  ~  2o(  <f>  —  <l>'  (68) 

02  —  =  0  +  0\  (69) 

and  for  the  D2R  term,  the  angles  are  the  same  as  for  RD2  with  0  and  <f>'  interchanged. 
The  RD  and  DR  fields  are  given  by 
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where  the  reflection  coefficient,  iZ«,;,,ts  ^  1. 

It  is  noted  that  the  term  D\Ris  present  in  the  region 
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h 

and  RD\  is  present  in  the  region 
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D-jR  is  present  in  tile  region 


arctan 


1-2  sin(a') 


[/t  —  I2  cos(a) 
and  RD2  is  present  in  the  region 

f  I2  sin(a) 
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The  case  where  a  is  90°  is  considered  first.  It  is  well  known  that  the  diffracted 
field  from  the  axis  of  the  dihedral  will  be  zero  for  this  case.  Using 

in  Equation  (70),  the  following  is  obtained 
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With  a  small  bistatic  angle  and  kl\  large,  the  second  term  is  dominant.  Thus,  the 
above  expression  may  be  written  as 

M  2/,  .  (<t>  +  <l>’\  {  si"  [2fcii  sin  (^)  sin  (^)] 
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For  the  backscatter  case,  with  strips  of  equal  length,  from  Equations  (72)  to  (75), 
only  D\R  and  RD\  are  present  for  the  range  0  <  ^  <  ir/4.  Also,  only  D2R  and  RD2 
are  present  for  the  range  7r/4  <  ^  <  ir/2.  For  the  backscatter  case 


and 


y"™’, -r/4  <  *  <  ,r/2. 


(79) 


(80) 
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Notice  that  for  <f>  =  45°  this  result  agrees  with  the  well  known  PO  result.  Figure  14 
shows  an  excellent  agreement  of  the  result  obtained  by  Equations  (79)  and  (80)  with 
that  obtained  by  the  UTD  computer  code  (with  all  terms  mentioned  above). 

Note  that  in  Figure  14  there  is  a  slope  discontinuity  in  the  dominant  UTD  (RD 


Figure  14:  The  dominant  UTD  terms  and  the  computer  code  results  for  backscatter 
for  9  in  90°  two-dimensional  dihedral  at  10  GHz. 


and  DR  terms)  at  ^  =  45°.  The  DD  field  across  the  face  of  the  dihedral  balances  this 
slope  discontinuity.  For  the  vertical  and  horizontal  polarizations,  Figures  15  and  16, 
respectively,  compare  the  UTD  with  and  without  the  DD  field  across  the  face  of  the 
dihedral  and  the  Method  of  Moments  solution  (18,  19]  for  the  backscatter  case  with 
if>  =  45°.  To  illustrate  this  point  further,  consider  the  4°  bistatic  case.  Figure  17 
shows  the  DR  and  RD  terms  for  a  4°  bistatic  angle.  Note  that  the  slope  discontinuity 
of  the  backscatter  case  now  spreads  out  and  makes  the  pattern  discontinuous.  This 
is  due  to  different  boundaries  of  the  various  RD  and  DR  terms  (Equations  (72)  to 
(75)).  Figure  18  shows  that  this  discontinuity  may  be  corrected  by  including  the 
DD  term  (see  Section  4)  across  the  face  of  the  dihedral. 
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DIAGONAL  LENGTH  IN  WL. 


Figure  15:  Contribution  of  the  DD  across  the  face  of  the  dihedral  for  the  backscatter 
case  at  45°  incident  for  the  vertical  polarization. 


DIAGONAL  LENGTH  IN  WL. 


Figure  16:  Contribution  of  the  DD  across  the  face  of  the  dihedral  for  the  backscatter 
case  at  45°  incident  for  the  horizontal  polarization. 
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Now  consider  the  case  when  n  is  not  close  to  00°.  It  is  expected  that  the  maximtim 
scattering  will  occur  when  4>~  <f)'  —  2a  —  — ir  or  (f>  —  (f>'  2a  —  w.  In  this  case,  either 
the  DR  or  RD  is  a  major  contributing  factor  along  with  the  diffraction  from  the  joint 
of  the  dihedral.  Notice  that  the  first  condition  is  satisfied  \i  <j>  >  <f>'  and  a  >  90°,  or 
<i>  <  <f>'  and  a  <  90°  {D\R  or  RD-j  term  is  the  major  contributing  factor  and  only 
one  term  is  present  in  a  given  region).  The  second  condition  is  satisfied  \{  <l>  >  <f>'  and 
a  <  90°,  or  <j>  <  <!>'  and  a  >  90°  (RDj  or  D^R  term  is  the  major  contributing  factor 
and  only  one  term  is  present  in  a  given  region).  Hence,  it  may  be  assumed  that  this 
is  a  region  where  diffraction  from  strip  1  and  reflection  from  strip  2  is  present.  It  is 
also  observed  that  when  the  DtR  term  is  large,  the  first  cotangent  term  of  the  singly 
diffracted  term  (Equation  (63))  is  also  large.  On  the  other  hand,  when  the  RDt  term 
is  large,  the  second  cotangent  term  of  the  singly  diffracted  term  (Equation  (63))  is 
large  too.  For  e  small,  the  case  when 


and 


~  2a  =  -  TT  +  e, 

( <f>  —  (/>'  —  2a\  e 

:os  ^  -  j  =  2’ 

’H — s — 


(81) 

(82) 

(83) 


is  considered.  First,  the  contribution  from  DxR  and  the  first  cotangent  term  of  the 
single  diffraction  may  be  simplified  to 


^ drmiinarti  m  i 


-2  jki, 
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Using  Equation  (81)  in  (85)  gives 
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Next,  the  rontribiitioTi  of  ftDx  and  the  second  cotangent  term  of  the  single  diffrac¬ 
tion  are  obtained  as 

_e-Wi 


j jin'll  _ 

2v/^ 


2i*/,  si„(la-l^i) 

sin  (2a  -I- 


i-(^) 


(87) 


At  this  stage,  to  simplify  Equation  (87)  further,  it  is  necessary  to  assume  that  a  is 
close  to  90°.  Consider  ^  to  be  a  small  angle  such  that 


a  =  7r/2  -f 

Now  Equation  (87)  may  be  written  as 

_  IT,  Jir/4  ) 

^(/ominnn/  ^ 

'4,^4,'-2a\  I  sin  [28kl,  cos 


(88) 
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\/A  V  2  J  \  2Skli  cos 
Again  using  Equation  (81)  in  (89)  gives 


r<i<  _  It,  -ji/ A  ~2jhl\tsU, {,!>') 
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(90) 


Had  the  i2Z?i  term  been  large,  the  result  would  have  been  the  same  as  Equations  (85) 
and  (89).  However  using  the  condition  for  maximum  scattering  for  this  case  {<f)  — 
<(>'  -\-2a  =  ir),  Equation  (85)  becomes 


jT0,U-\-0\ft,  f.„,  fji  in/ 1  •  fr,  i,\ 

Ud,„„i„ani  =  ^  Sin  (2a  -  0  ) 


(91) 


and  (89)  becomes 
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^  ^  .m,h^  ^ 
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(93) 

(94) 

(95) 


wjf^r>‘i  +  0‘2nH,eat  _  Tp  — 27^/56  Rin(^'—a) 

^Hominanl  ^ 


sin  [26kl2  sin  {(f/  —  a)] 
2Skl2  sin  {<f/  —  a) 


!■ 


(96) 


It  is  important  to  emphasize  that  Equations  (84)  to  (96)  are  present  over  regions 
as  given  by  Equations  (72)  to  (75)  and  are  valid  for  one  of  the  two  cases  of  max¬ 
imum  scattering  as  mentioned  earlier.  Figures  19,  20  and  21  are  for  equal  length 
strip  dihedrals  of  a  =  92°  (4°  bistatic),  95°  (10°  bistatic)  and  110°  (40°  bistatic),  re¬ 
spectively,  and  compare  the  results  obtuned  by  the  dominant  UTD  terms  with  that 
obtained  by  the  UTD  computer  code.  Good  agreements  are  obtained  in  all  regions 
except  where  the  speculars  from  the  strips  are  present,  which  are  not  included  in 
the  dominant  UTD  terms. 


3  Specular  Return  From  Dihedral  Strips 

As  pointed  out  earlier,  the  specular  from  the  individual  strip  faces  is  produced  by 
D  from  the  ends  of  the  strips  and  the  DR  and  DRD  terms  between  the  strips. 

The  importance  of  D  and  the  DR  terms  is  illustrated  by  the  following.  Consider 
the  geometry  shown  in  Figure  22.  It  is  observed  that  the  singly  diffracted  field 


39 


Figure  19:  The  dominant  UTD  terms  and  the  UTD  computer  code  results  for  9  in 
92°  2-D  dihedral  at  10  GHz  (4°  bistatic). 


0  (deg.) 


Figure  20:  The  dominant  UTD  and  the  UTD  computer  code  results  for  9  in  95°  2-D 
dihedral  at  10  GHz  (10°  bistatic). 
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Figure  21:  The  dominant  UTD  terms  and  the  UTD  computer  code  results  for  9  in 
110°  2-D  dihedral  at  10  GHz  (40°  bistatic). 


R  D, 


Figure  22:  Cancellation  of  infinities  of  diffracted  and  diffracted-reflected  terms. 

D\  appears  to  leave  an  infinity  at  the  reflection  shadow  boundary.  The  diffracted 
-  reflected  field  D2R  also  appears  to  leave  an  infinity  along  the  incident  shadow 
boundary.  However  and  D2R  combine  to  give  a  bounded  result. 

The  diffracted  field  is 


(97) 


where 


and 


Kh  =  < 


E',  £  parallel  to  wedge 

H',  E  perpendicular  to  wedge 


(98) 


/  1  1  \ 

icos  j  ’ 

where  the  reflection  coefficient,  /?,,/,  =  ^1.  Along  the  reflection  shadow  boundary, 

=  cos(ir/2  -  e/2)  =  e/2.  (100) 


cos 
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For  the  DR  field, 


(101) 


U:,,R,J,D,J.(Q.2)e 
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where 


Dn.hiQ'i)  = 


/  \  "t"  " 


2v^  \^cos{^)  ’  ‘cos  (^)  /  ■ 


(102) 


Along  the  incident  shadow  boundary, 

cos  ^^  '2  “  cos(ir/2  +  e/2)  =  —e/2.  (103) 

By  adding  the  singular  terms  of  D\  and  D2R,  the  following  is  obtained 


-e--”'/'  /-2 


where  the  phase  contributions  (see  Figure  23)  are 


(^) 


,  (104) 
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and 

=  e-.l*rg-jM./+(^+0ros(«+r))_ 


(106) 


Equation  (104)  can  be  written  as 
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(107) 


As  e  approaches  zero.  Equation  (107)  becomes 

fj„nnulnr  ,n,nt  ^  {kl3in{d>)}  .  (108) 

\/2nk 

Note  that  this  gives  the  PO  result  for  a  strip  of  length  1. 

It  is  observed  that  for  the  backscatter  case  from  a  90°  dihedral,  the  DR  term  is 
shadowed  due  to  the  diffraction  reflection  shadow  boundary  —  90°).  As  pointed 
out  in  the  next  section,  it  is  the  DRU  term  that  takes  care  of  this  discontinuity. 
Figure  24  shows  that  the  inclusion  of  DRD  term  ensures  the  continuity  and  gives 
the  correct  result.  Note  that  al  (^  =  90°  only  the  9  in  strip  is  visible,  which  by  PO 
gives  a  result  of  10.39  dBm.  The  result  obt lined  by  Figure  24  is  10.26  dBm. 
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Figure  24:  Importance  of  DRD  term.  The  results  are  for  a  9  in  90°  2'D  dihedral  at 
10  GHz. 


4  Continuity  of  Fields  by  Including  Double  Diffrac¬ 
tion 

Figure  25  shows  the  shadow  boundaries  that  can  result  after  the  first  diffraction. 
The  diffracted  field  is  discontinuous  across  the  diffraction  incident  shadow  boundary 
(DSB),  and  the  diffracted  -reflected  field  is  discontinuous  across  the  diffraction  reflec¬ 
tion  shadow  boundary  (DRSB).  The  DxD2  field  compensates  for  the  discontinuity  of 
D\  and  D\  R  fields.  The  double  diffraction  formulation,  for  separate  wedges,  as  given 
by  Equation  (51),  has  16  terms.  Across  the  DSB  or  DRSB  one  of  the  terms  becomes 
infinite.  This  term  is  determined  by  the  values  of  $p  and  (Equations  (49)  and 
(50))  which  make  the  cotangent  (cot  and  cot  (^))  unbounded.  Thus,  any 
term  can  become  singular  depending  on  the  geometry  and  the  position  of  the  source 
and  the  receiver.  As  shown  in  an  example  below,  the  same  holds  true  if  the  wedges 
are  joined. 

Consider  the  case  when  the  two  wedges  are  not  joined,  as  shown  in  Figure  26. 
Since  ir  —  (^i  —  ^\)  =  0,  cot(O)  =  oo,  and  ir  {(f>2  -f  ^2)  =  4n,  cot(7r)  =  00,  at  the 
boundaries,  thus  the  term  P($,,=  i,#,=4)  will  balance  the  singularity  of  the  singly 
diffracted  field.  The  terms  F(#,,r=i ,4>,=  i)  -f  F(^p=i,$,=2)  +  /’($,,=  ], $,_t)  balance 
the  discontinuity  that  results  from  the  diffraction  by  the  edge  Q2  being  shadowed 
by  edge  Qi.  The  terms  F($p=2,  ^,=4)  +  F(4*,,=,i,  ^,=j)  +  F($p=, ,$,=  ,)  balance  the 
discontinuity  that  results  from  the  diffraction  by  edge  Qt  being  shadowed  by  edge 
<?2- 

In  the  above  sections,  the  doubly  diffracted  field  was  necessary  to  get  the  con¬ 
tinuity  of  the  fields  in  two  cases.  (See  Figures  24  and  18.)  In  both  of  these  cases 
the  discontinuity  is  a  result  of  being  on  the  DRSB.  In  this  section,  a  case  where  the 
discontinuity  is  along  the  DSB  is  considered.  It  is  shown  mathematically,  as  well 
as  by  plots,  how  the  inclusion  of  doubly  diffracted  field  compensates  for  the  singu¬ 
larities  and  discontinuities.  In  Figure  27,  C)  ~  0  and  62  ~  0  i.e.,  both  incidence 
and  observation  aspects  are  close  to  grazing.  Note  that  for  ci  <  0  and  £2  >  0, 


46 


Figure  25:  Illustration  of  diffraction  incident  (DSB)  and  diffraction-reflection  shadow 
boundary  (DRSB). 


Figure  26:  Balancing  of  fiehls  for  disjoined  wedges. 


diffrartion  from  D\iDi  and  is  present  But  for  f|  <  0  and  tj  <  0,  only 

diffraction  from  D-i  and  is  observed.  Details  of  the  equations  below  are  given  in 
Appendix  A.  The  contribution  of  the  fields  singly  diffracted  by  Qi  can  be  expressed 
as  a  singular  contribution  given  by 


■\/2irk 


e  [1  +  «sn(e2)] 


and  a  non-singular  contribution  of 


(109) 


(110) 


Similarly,  separating  Pe,miQ2)  into  a  singular  and  a  non-singular  part,  where  the 
singular  part  is  given  as 


KM)  = 

and  the  non-singular  part  is 
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For  the  single  diffracted  field  from  Q;»,  the  singular  term  is  given  as 

-e-J’'/'  r  1 


K.M)  = 


\/2nk 


jw/i  1 
Trfc  Ici  —  C2 


(113) 


The  double  diffraction  term  obtained  by  setting  p  =  1  and  g  =  1  in  Equation  (44) 
compensates  for  the  discontinuity  and  singularities  of  the  single  diffracted  fields.  This 
term  is  given  as 

.g-W  f  ev5gn(e2)  -  e,  5511(6, )'» 


P t,h  (  I  >  ^7=  1  )  ~ 


y/2'Kk 


6T  -  6; 


1  ±  1 
- e 


-  jkti 


■)kfi 


(114) 
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Adding  all  the  singular  terms  of  singly  diffracted  fields  from  Qi,  Q2  and  Q:,,  and 
doubly  diffracted  field  from  Qi  and  Q2,  Equations  (109),  (111),  (113)  and  (114),  the 
following  is  obtained 


n.^fffgular  total  ^  ^ _ 

\/2^  k 


gjkh(ti  -<j)' 

Cl  —  €2 


1  ±  1  \/d 

2  IT  ' 


(115) 


Now  for  6 1  and  62  zero,  the  above  equation  becomes 


pningutar  total  _  _ 


\±ly/d  _ 


jkd 


(116) 


Note  that  the  first  term  is  equivalent  to  a  PO  result  of  the  strip  of  length  h. 
The  second  term  is  the  correction  to  the  PO.  Figure  28  shows  that  the  inclusion 
of  the  diffracted-diffracted  term  ensures  continuity.  Thus,  the  UTD  formulation  of 
the  double  diffraction  presented  in  (2)  provides  the  compensation  of  the  singulari¬ 
ties,  so  that  the  leading  term  of  the  total  scattered  far  field  is  well  behaved  in  the 
neighborhood  of  the  shadow  boundary  (DSB)  of  the  diffracted  incident  field. 

Now,  consider  the  continuity  of  the  non-singular  part  of  the  diffracted  terms.  It 
is  also  noted  that  vanishes  due  to  shadowing  by  the  edge  Q2  (c2  <  0), 

and  similarly,  ■P"„,(<?2)  vanishes  due  to  shadowing  by  the  edge  Q\  (c|  >  0). 

From  the  discussion  in  the  beginning  of  this  section,  it  is  expected  that  the  term 
Ph[^p=-2i^q=\  —  C2)  ensures  continuity  of  Pr[m(Q\)-  From  Equation  (110),  when 

and  from  Equation  (42)  and  (43), 


^p=2  =  ’T  -f-  (n,7r  -  (n|7r  -  TT  -  c,))  =  27r  -f  e,  (118) 

=  n  -  (n-jir  ~  (tizTT  -  tt  -  cj))  =  ^2  (119) 

N,=2  =  =  1  for  0.67  <  n  <  2  (120) 

2n7r 

Op-j  =  2cos^  [(27r  -I-  e|)/2  -  n/2  -  2n7-/2]  —  2sin‘  [7r(l  -  n)  -t-  ei/2]  (121) 
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Figure  28:  Importance  of  the  doubly  diffracted  term  for  the  field  to  be  continuous 
for  TE  case. 
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again  neglecting  terms  of  order  cf  and  ei|, 

(^) 

Thus  P,„  +  Pf,  (Equations  (117)  and  (123))  is  continuous  across  the  shadow  bound¬ 
ary.  Similarly  it  may  be  shown  that  the  term  P/i(^p=i,^,r:2)  ensures  continuity  of 

KiQ^y 

However,  for  a  TM  plane  wave  illumination,  P^'(Qi)  and  P/.'{Q2)  vanishes  rapidly 
at  C)  =  C2  =  0,  so  that  no  discontinuity  occurs  due  to  the  shadowing  by  Qi  and  Q2, 
However,  the  terms  /*,($,,=2j^q=i)  and  P,($p=i,$,=2)  provide  a  contribution  which, 
as  expected,  is  continuous  across  the  DSB  and  provides  a  smoother  transition  from 
the  lit  into  the  shadow  region  of  the  field  for  rays  singly  diffracted  from  edges. 

From  Equation  (45), 


P n{<f>p=2td>q=\)  = 


2ell4 


87ry(n,n2)q,i„.  (^)  (^) 


(124) 


=  L-f _ \ _ 

SirynJ 

Figure  29  shows  how  the  result  is  affected  by  including  the  DD  field  for  the  case 
when  the  TM  plane  wave  is  incident. 


5  More  Multiple  Plate  UTD  Terms 

The  previous  sections  discussed  the  plate  interaction  terms  that  provide  the  ba¬ 
sic  dihedral  characteristic  pattern.  This  section  briefly  outlines  other  higher  order 
mechanisms  that  contribute  to  the  pattern.  Depending  on  the  angle  of  the  dihedral 
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Figure  29:  For  the  TM  case,  the  doubly  diffracted  field  provides  a  smoother  transi¬ 
tion. 


and  on  the  bistatic  angle  any  one  of  these  terms  ran  be  important.  In  general,  ns 
the  frequency  decreases,  that  is,  as  the  size  of  the  object  in  wavelengths  gets  smaller, 
more  higher  order  terms  need  to  be  added. 

The  various  UTD  terms  can  be  derived  by  using  combinations  of  the  theory  in 
Section  2.  The  importance  of  single  diffraction  (D),  DR  and  RD  terms  in  giving 
rise  to  the  main  beam  of  the  pattern,  and  the  single  diffracted,  DR  and  the  DRD 
terms  in  giving  the  specular  from  the  strips  has  already  been  considered.  The  RRD 
and  the  DRR  terms  are  present  when  the  angle  of  the  dihedral  is  less  than  90°.  Of 
course,  these  two  terms  are  most  significant  when  the  diffraction  is  along  a  reflection 
shadow  boundary  RSB.  Although  the  RDR  term  is  present  for  all  angles  of  the 
dihedral,  its  significance  is  greatest  when  the  angle  of  the  dihedral  is  less  than  90° 
and  the  diffraction  is  along  a  reflection  shadow  boundary. 

A  discontinuity  will  result  if  a  reflection  falls  off  the  edge  for  which  a  diffraction 
must  be  included  to  balance  the  fields.  Therefore,  the  DDR  term  compensates  for 
the  discontinuity  of  the  DRR  and  RDR  terms.  The  RDD  term  compensates  for 
the  discontinuity  of  the  RRD  and  RDR  terms.  Also  the  DRD  compensates  for  the 
discontinuity  of  the  RRD  and  DRR  terms.  Discontinuities  of  the  DDR,  DRD,  and 
RDD  are  balanced  by  the  triply  diffracted  fields.  Thus,  one  term  may  be  balanced 
by  a  number  of  different  terms  and  vice  versa  depending  on  the  geometry. 

6  Results  of  Rectangular  Dihedral 

In  this  section,  the  theory  and  computer  code  are  validated  by  comparing  the  UTD 
results  to  the  results  based  on  a  2-D  method  of  moments  [18,  19]  solution  of  the 
rigorous  integral  equation  for  the  currents  induced  on  the  target  by  the  illumination. 
Unless  explicitly  shown  in  the  pattern,  the  geometry  of  the  dihedral  is  as  shown  in 
Figure  12.  In  order  to  calculate  the  backscatter  results,  a  slight  bistatic  angle  is 
used  to  avoid  numerical  limit  difficulties  w'hen  approaching  indeterminate  forms 
on  the  computer;  alternatively,  one  could  calculate  the  limits  analytically  using 
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L’HopitaPs  rule  and  calculate  the  fields  without  requiring  one  to  use  slight  lusiafic 
angles.  Method  of  Moments  is  chosen  to  validate  the  results  as  an  independent 
check,  since  both  are  two-dimensional.  Fairly  large  sized  dihedrals  can  be  calculated 
using  Method  of  Moments  in  the  2-D  case.  In  3-D,  limits  on  computer  memory 
restrict  the  size  of  dihedral  that  can  be  analyzed  using  Method  of  Moments.  A  9  in 
2-D  dihedral  is  analyzed  at  2  GHz  and  10  GHz  which  means  the  dihedral  is  1.52 
and  7.62  wavelengths  on  a  side,  respectively.  All  the  patterns  are  normalized  to  dB 
below  a  meter,  that  is  a  two  dimensional  cross  section. 

For  the  vertically  polarized  incident  field,  the  backscatter  from  a  90°  dihedral. 
Figure  30  shows  the  singly  diffracted  (D)  term,  DR,  and  RD  terms.  Figure  31  shows 
the  D,  RD  and  DR  terms  combined,  the  DRD  term,  and  the  DD  term.  This  is  all 
that  is  required  to  get  the  complete  result  for  the  90°  dihedral.  Figure  32  shows  the 
result  obtained  by  the  D,  DR,  RD,  DRD,  and  DD  terms  combined.  Notice  that  the 
result  obtained  by  Figure  32  is  the  same  as  that  obtained  by  all  the  terms  of  the 
UTD  computer  code.  (See  Figure  45.) 

For  the  backscatter  case  with  strips  of  equal  lengths,  based  on  simple  reflection 
concepts,  the  order  of  interactions  required  to  get  the  pattern  in  between  the  strips  is 
Where  INT  denotes  the  integral  value  of  Thus,  for  q  =  60°  third 

order  terms  are  required  to  get  the  pattern  in  between  the  strips.  Figure  33  shows 
how  only  the  third  order  terms  for  a  dihedral  of  60°  combine  to  give  the  correct  result 
in  between  the  strips.  Figure  34  compares  the  Method  of  Moments  results  with  that 
of  the  UTD  computer  code  (with  all  terms  present)  for  a  dihedral  of  60°.  Notice  that 
the  discontinuity  at  0  =  60°  in  Figure  33  and  Figure  34  is  due  the  triply  diffracted 
field  which  has  not  been  included  in  the  computer  code.  Our  formulation  indicates 
that  for  a  =  45°,  fourth  order  interactions  are  required.  Since  the  UTD  computer 
code  does  not  take  into  account  all  the  fourth  order  interactions,  we  expect  no 
agreement  of  the  results  in  between  the  strips.  This  is  shown  in  Figure  36.  However, 
even  for  this  case,  notice  the  excellent  agreement  for  (f)  greater  rhen  50°.  Thus,  for  the 
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Figure  30:  D,  DR,  and  RD  terms  for  the  backscatter  for  9  in  90°  rectangular  dihedral 
at  10  GHz  (vertically  polarized). 
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Figure  31:  All  terms  of  the  above  figure,  DRD,  and  DD  terms  for  the  backscatter 
for  9  in  90°  rectangular  dihedral  at  10  GHz  (vertically  polarized). 
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Figure  3^  All  terms  required  to  get  the  backscatter  result  for  9  in  90°  rectangular 
dihedral  at  10  GHz  (vertieally  polarized). 


terms  present  in  the  rompnter  rode,  the  pattern  would  worsen  graref»ll_v  in  between 
the  strips  for  a  less  then  60°.  As  the  angle  of  the  dihedral  becomes  smaller,  more 
terms  are  required  to  give  the  correct  result  in  between  the  strips  and  a  “cavity” 
effect  results. 

From  here  on,  the  UTD  results  contain  all  of  the  terms  discussed  above.  For  the 
backscatter  case,  from  Geometrical  Optics,  a  main  beam  is  expected  when  the  angle 
of  the  dihedral  is  such  that  is  an  integer.  This  is  clearly  evident  in  Figures 

32,  35  and  36.  Also,  as  the  angle  of  the  dihedral  diverts  from  a,  which  makes 
an  integer,  the  main  beam  falls.  Figures  37  to  41  show  backscatter  results  for  a  = 
85,  92,  94,  95,  and  120°  with  the  incident  field  horizontally  polarized.  As  expected, 
the  main  beam  of  the  pattern  for  a  =  90°  falls  as  a  is  decreased  or  increased. 

The  backscatter  field  for  the  2  GHz  case  is  shown  in  Figure  42  for  horizontal 
polarization  and  in  Figure  43  for  vertical  polarization.  The  10  GHz  backscatter 
case  is  shown  in  Figures  44  and  45  for  the  horizontal  and  vertical  polarizations, 
respectively.  Patterns  for  a  fixed  bistatic  angle  of  90°  are  shown  next.  The  2  GHz 
case  is  shown  in  Figure  46  for  horizontal  polarization  and  in  Figur#^7  for  vertical 
polarization.  The  10  GHz  patterns  are  shown  in  Figures  48  and  49  for  horizontal 
and  vertical  polarizations,  respectively.  In  all  cases  the  agreement  is  excellent.  The 
small  dihedrals  have  slightly  more  disagreement  due  to  the  fact  that  the  higher 
order  interactions  are  more  important.  It  is  assumed  that  even  better  agreement 
and/or  smaller  sized  dihedrals  could  be  achieved  by  including  more  interactions. 
For  engineering  accuracy,  however,  this  would  not  seem  to  be  necessary  in  this 
application. 
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Figure  33:  Only  third  order  terms  give  continuous  result  for  9  in  60°  rectangular 
dihedral  at  10  GHz  for  horizontally  polarized  field  (backscatter). 
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Figure  34:  Backscatter  field  from  a  9  in  60  rectangular  dihedral  at  10  GHz  for 
horizontally  polarized  field. 
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Figure  35:  Backscatler  result  tor  9  in  60°  rectangular  dihedral  at  10  GIU  for  hon- 
zontally  polarized  field. 


Figure  36:  Backscatter  field  from  a  9  in  45”  rectangular  dihedral  at  10  GHz  for 
horizontally  po  arized  field. 
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Figure  39:  Backscatter  result  for  9  in  94°  rectangular  dihedral  at  10  GHz  for  hori¬ 
zontally  polarized  field. 


Figure  40:  Backscatter  result  for  ■)  in  95"  rectangular  dihedral  at  10  GHz  for  hori 
zontally  polarized  field. 


Figure  41:  Backscatter  result  for  9  in  120°  rectangular  dihedral  at  10  GHz  for 
horizontally  polarized  field. 


Figure  42:  Backscatter  field  from  a  9  in  rectangular  diliedral  at  2  GHz  for  horizontally 
polarized  field 
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Figure  43:  Backscatter  field  from  a  9  in  rectangular  dihedral  at  2  GHz  for  vertically 
polarized  field. 


Figure  44:  Backscatter  field  from  a  9  in  rectangular  dihedral  at  10  GHz  for  horizon¬ 
tally  polarized  field. 
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Figure  49:  Ninety-degree  bistatic  field  from  a  9  in  rectangular  dihedral  at  10  GHz 
for  vertically  polarized  field. 
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SECTION  4 


Rim  and  Curved-Surface  Caustic 
Corrections 

1  Introduction 

The  Geometrical  Theory  of  Diffraction  predicts  infinite  fields  at  ray  caustics.  The 
occurrence  of  a  caustic  requires  that  the  erroneous  ray  field  infinities  be  “corrected” 
with  transition  functions  which  blend  the  well  behaved  wide-angled  ray  field  predic¬ 
tions  smoothly  into  the  corrected  or  bounded  values  at  the  caustic.  Details  of  rim 
and  curved-surface  caustic  corrections  can  be  found  in  References  [20,  21,  22].  Al¬ 
though  theses  solutions  are  based  on  equivalent  edge  and  line  currents,  the  modern 
MEC  (Michaeli  et  al)  provides  essentially  the  same  result. 

Ryan  and  Peters  [20]  developed  a  solution  in  the  rim-caustic  regions  based  on 
the  field  radiated  by  equivalent  edge  currents  on  the  rim.  This  result  goes  to  the 
diffraction  point  solution  outside  the  caustic  region;  thus,  it  is  a  uniform  solution. 
In  Reference  [21],  this  technique  is  extended  by  the  use  of  the  “stripping  concept” 
[23,  24]  in  addition  to  using  a  small  argument  approximation  to  the  diffraction 
coefficient  in  the  specular  caustic  region  of  the  radiation  integral.  This  solution 
referred  to  as  the  “Bessel-Struve  function  extension,”  provides  an  accurate  uniform 
result  in  the  rim-caustic  region  that  is  associated  with  a  specular  or  forward-scatter 
direction  of  disks,  end  caps  of  cylinders,  cones  etc. 
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The  need  for  the  rnrved-surfare  specular  correction  arises  from  the  fact  tltaf 
the  classical  edge  diffraction  coefficient,  being  a  local  phenomenon,  has  no  a  priori 
knowledge  of  the  variation  in  curvature  along  a  curved-surface  away  from  the  edge. 
If  the  radius  of  curvature  is  equal  at  both  the  edges  of  the  curved-surface,  as  in 
the  case  of  the  cylinder,  the  total  field  along  the  specular  direction  will  be  correct. 
However,  if  the  radii  of  curvature  are  different  at  the  two  edges  of  the  curved-surface 
(for  example,  cones),  the  energy  spreads  differently  at  every  point  on  the  curved- 
surface,  depending  on  the  radii  of  curvature.  Since  this  information  is  not  built 
into  the  classic  edge  diffraction  coefficient,  the  edge  diffracted  field  is  not  able  to 
match  into  the  reflected  field  in  the  specular  direction  whose  dominant  behavior  is 
as  predicted  by  physical  optics  (PO).  Using  procedures  similar  to  References  [20,  21) 
and  extending  them  by  asymptotically  matching  the  PO  solution  for  the  reflection 
from  the  curved  face  into  the  UTD  solution,  Ebihara  [22]  extracted  a  transition 
function  useful  in  the  curved-surface  specular  region. 

In  Section  4.2,  the  equivalent-line  current  solution  procedure  [25]  is  presented. 
In  Section  3.1,  the  equivalent-edge-current  solution  of  Ryan  and  Peters  is  presented. 
Next,  in  Section  3.2,  the  Bessel-Struve  function  extension  is  examined.  Finally,  the 
curved-surface  specular  correction  will  be  considered  in  Section  4. 

2  Equivalent-Line  Current  Solution  Procedure 

The  following  steps  are  used  to  find  the  equivalent-line  solution: 

1.  The  field  of  a  canonical  problem  is  obtained  (a  two-dimensional  solution). 

2.  Next,  the  radiated  fields  from  infinite  travelling-wave  line  currents,  I''  and 
are  determined.  I'  is  an  electric  current  and  is  a  magnetic  current. 

3.  The  currents  /'  and  V”  are  substituted  for  the  two-dimensional  solution, 
and  the  currents  I'  and  /"'  are  then  adjusted  to  give  f  lie  same  field  as  the  two- 
dimensional  solution.  In  the  principal  plane,  /'  gives  the  horizontally  polarized 
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fieUl,  anil  P"  gives  tlie  vertirally  imlariveil  firlil.  Tliese  are  “firlif.ious”  rnrrenf.s  that 
vary  with  the  incident  and  observation  positions. 

4.  These  equivalent  line  currents  are  broken  up  into  infinitesimal  current  ele¬ 
ments.  The  current  elements  are  placed  at  the  equivalent  locations  on  the  surface. 

5.  The  radiated  fields  due  to  the  above  currents  are  evaluated.  This  involves 
a  line  integral.  The  resulting  expression  is  the  equivalent-line  current  solution.  It 
gives  the  correct  value  of  the  field  away  from  the  caustic  direction. 


3  Rim  Caustic  Correction 

3.1  Ryan/Peters  Equivalent-Edge- Current  Solution 


Ebihara  [22]  used  the  original  work  by  Ryan  and  Peters  [20]  to  get  the  standard 
diffraction  coefficient  with  rim  caustic  correction  function,  Figure  50.  There  appear 
to  be  some  sign  discrepancies  early  in  the  derivation  in  [22]  which  result  in  incon¬ 
sistencies  in  the  final  result  for  the  hard  polarization  using  Equation  5.225  in  [22]. 
In  this  section,  the  standard  diffraction  coefficient  with  the  transition  functions  are 
defined  for  consistency. 

The  equivalent  edge  current  used  by  Ryan  and  Peters  are 


and 


where 


and 


r  = 


2jG‘ 


T}k  sin  /3„  sin 


7  (f  ■  ^'(0))  ' 


jkf,  r' 


2jG" 


k  sin  (3„  sin  /?' 


(c.ff'(O)) 


lArr,  .r' 


C”  (*,♦•)  =  4  B) 


.4  (#,  #')  =  cot 


/tt  -f  ($  - 
\  ^ 


2n  }  ' 


(126) 

(127) 

(128) 

(129) 

(130) 
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JB  =  cot  ^ - - - -j  +  cot  ^ - 


t  'f*' 


(131) 


These  currents  are  distributed  on  the  illuminated  edge  of  the  rim.  The  radiated 
fields  due  to  the  electric  and  magnetic  equivalent  edge  currents  on  a  circular  rim  are 


a  cos  6  e  \  t-  l 

- ~L  (132) 

an  f  G’"  /  -  v  ,, 

^ -  /  -^(e..ff'(0))cos(.^-<A,.)e^^"-(^+^)d.^.,  (133) 

lir  r  Jit  sin  p„  '  ' 

=  (e-E’(0))cos(<^-.^,)e^"^"'  (^+^')d^^^^  (134) 

ZTT  r  Jtt  sin  47„  '  ' 


ET  = 


27r  r 


^  (e-^'  (0))sin(,^-,^,,)e^'^^''-('+^'W..  (135) 

27r  r  Jit  sin  p„  ^  ' 

The  diffracted  field  is  the  sum  of  the  fields  radiated  from  both  electric  and  magnetic 
equivalent  currents,  so  that 

Ei„,  =  0Eo  +  ^E^,  (136) 


where 


E^t,  =  £:;  +  e: 


(137) 


E„  =  E’o  +  e;;’. 


(138) 


To  simplify  the  integration  in  (132)-(135),  Ryan  and  Peters  used  the  approxima¬ 
tion  that  the  diffraction  coefficients  remained  relatively  constant  with  respect  to  the 
variable  of  integration,  so  that 


sin  /3,. 


constant  in  (f), . 


(139) 


V»,  -  kah,  (</>, )  •  (r  4  >  ') 

Icr  [sin  0\  os  (</>’  c/>, )  4  .sinfl  cus  (</>  -  </.  )| 


(140) 

(141) 


and 


7.  =  kzj  •  (f  +  f')  . 

Using  the  previously  mentioned  approximation  and  the  fact  that 

/  sin  <f),  cos  <f>re^'^' d<f>r  —  0, 

Jo 

the  expressions  for  the  radiated  fields  become 


a  e-’’*' 

— - E’Jcr  co^^4>.e^'‘'d4>, 

ZTT  Sin  p„  r  I  Jn 

—  G'”  cosffcosO'J^^  sin^  d(f)f^  , 


E„  =  — 


CL  C 

- — : — 2 - E'o  ^  G"'  COS  8  cos  8'  I  sin^  ■f  t  '  ■  IS, 

2ir  sin  Po  r  {  Jt) 

-G"'  'dr,  j, 


where 


Vv  -  U  cos  <!>,■. 


Let  the  following  integrals  be  defined  by: 


Fi{V)—  f  cos' (p,  e^'" d<f>, 

Jn 

/,((/)  =  J'  cOS^d>re^'^'dd>r 

/,*((/)  =  I  '  cos^  ’ d<i>, 

Fi  {U)  —  cos  6  cos  6'  f  sin*  (f>,  c^'"  dd), 
Ji> 

f-2{U)  —  ^os  6  cos  8'  J  s\n' d>,  e^'^  ' d<j>, 

V 

fj  {U)  —  r  ns  8  cos  8'  j  sin'  (j),  f '  dd>, 
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(148) 

(149) 
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These  inlegrnis  ran  he  evalnafed  in  closed  form  as  follows 
/. 


(t/)  =  »  |(j,.((/)  -  +}  (»,({/)  -  (153) 

/;([/)  =  :r  |(y„(f/)  -  ■h^'j  -  j  (h„(U)  -  (154) 


/,(£!)  = 


(l-)5) 


/2'(U)  =  ircosScosfl' ,  (156) 

where  J|  and  Hi  are  the  first  order  Bessel  and  Hankel  functions,  respectively.  The 
expressions  for  the  radiated  fields  are 
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a  e^'"  e~^^’  1 

E},—  {~{A-B)F,  +  {a  +  B)Fi}. 


27rsin/3„  r  2n 
The  above  can  be  rearranged  into 

Breaking  the  integral  apart  into 

E,{u)  =  fi{u)  +  f;{U) 

F2{U)  =  h{u)  +  f;{u) 

and  further  manipulating  them  gives 
-a  e^'>'  I 
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Now,  consider  associating  a  half-rim  integral  with  the  points  at  <!>,  —  0,  tt,  the 
stationary  points  for  a  large  they  are  also  the  diffraction  points  on  the  rim  for 
the  principal  plane.  From  the  above  equations,  the  field  contributions  from  the  half 
rim  associated  with  <f>c  —  0  are 

—ft 

-  ^(/'  +  fv}  (165) 

27r  sinpo  T" 

E^o  =  {A{f,  -  h)  +  B  if,  +  f,)}  (166) 

2a-  sin  Po  T" 

^^E,^  +  eEi„.  (167) 

Comparing  this  with  the  expression  for  the  edge-diffracted  field  from  a  point  on  the 
rim, 

Ef.  =  £’(Q.)\/p''(Qr)e-"^-"~,  (168) 

the  new  diffraction  coefficients  may  be  expressed  as: 


where 


where 


and 


D{Q,)  =  -0'm^iQr)  - 
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2nv27rfcsmp„ 

(170) 

2n\/2nk  sin  p„ 

(171) 

T..m  =  i  (/,((/)  -  /,((/)) 

(172) 

TMV)  =  1  (Ml/)  +  Ml/)) 

(173) 

[/  (Q,  )  ~  ka,n,  ■  (r  +  f'). 

(174) 

The  functions  T,.i  and  T,2  are  the  transition  functions  through  the  rim  caustic  region 
based  on  the  Ryan/Peters  equivalenl-edge-current  solution.  The  transition  functions 
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have  small  and  large  argument  fr>rms  whieh  are  obtained  using  t)>e  series  and  the 
asymptotic  forms  of  the  Bessel  and  Struve  functions,  so  that  for 
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Associating  a  half  rim  to  the  stationary  points  Q\  and  Q2,  the  total  field  radiated 
by  the  equivalent  edge  currents  on  the  rim  may  be  written  in  the  diffraction-point 
form 


F'’  F'’  I  F'’ 


(179) 
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where 


E'l  =  — 

w  ■p 

and 

Ei  =  E' 

T 

3.2  Bessel-Struve  Function  Extension 


(180) 


(181) 


One  problem  associated  with  the  Ryan/Peters  result  is  that  the  approximation  (139) 
becomes  invalid  for  a  rim  caustic  associated  with  an  end  cap’s  specular  or  forward- 
scatter  direction.  To  overcome  this  difficulty,  a  modified  edge  vector  e*  is  used.  It 
is  defined  by 

e'=p(ep),  (182) 


where 


P  = 


(f'  +  r)>i  n„ 


(183) 


j(f'  4-  f)  X  n„) 

n„  =  normal  to  the  end  cap 
f'  =  direction  of  source,  and 
r  ~  direction  of  receiver. 

By  using  the  modified  edge  vector,  only  the  perpendicular  components  to  the 
plane  of  incidence  of  the  equivalent  edge  currents  are  considered.  This  is  called  the 
“stripping  concept.”  The  equivalent  edge  currents  with  this  modification  are  given 


by 


K  = 

IT  = 
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Tfk  sin  sin  /?' 
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(184) 


(185) 


fe  sin  0„  sin 

The  modified  equivalent-edge-current  solution  is  obtained  in  the  same  way  as  before. 
To  simplify  the  evaluation  of  the  integral,  the  following  small  argument  approxima¬ 
tions  are  used: 

0  ^  in  the  specular  region^  or  (186) 
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6  ^  T  —  0'  —  e  ,  in  the  forward  scatter  region. 


(187) 


This  results  in 


and  the  integral 


Q,,.n  1 
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sin  (3„  sin  /?'  c  cos  (f>^ 
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which  can  be  evaluated  by  using 


f4U)  =  pCOB<f>.e^^’^"^'^^d4>r 


(188) 
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Associating  a  half  rim  to  the  stationary  points  Q\  and  Q2t  the  total  held  radiated 
by  the  equivalent  edge  currents  on  the  rim  may  be  written  in  the  diffraction-point 
form 

E;',„  =  Ef+£?  (191) 

where 

Ei  =  (192) 

Ei  =  £'J(C?2)v//(Q2)e-’*'^--'— .  (193) 

r 

The  modified  diffraction  coefficient  is 


D{Qr)  =  -4W(Qc)  -  4'^A.((?c)  (194) 


where 


^»,/i  = 


_#,-W 


2nv/SS 


P  =  =  d>±<t>\ 


(195) 


(196) 


V  =  kah,  •  (f  f') , 


(197) 


78 


(198) 


and 

fo{U)  =  i^\^- H^iU)  -\^jJ^iU)Y  (199) 

J\  and  H\  are  the  first  order  Bessel  and  Struve  functions,  respectively.  The  transition 
function  is  T,  (U),  k  is  the  wave  number,  a  and  are  the  radius  of  curvature  and 
the  normal  to  the  edge  at  the  point  of  diffraction,  and  f  and  f'  are  the  directions 
to  the  source  and  the  receiver,  respectively.  The  transition  function  has  small  and 
large  argument  forms  given  as: 
for  |[/|  <<  1 


for  \U\  »  1 

T  (t/)  ^  1  +  _  l./Ae-.>("-7)  +  _i^ _ il--.  (201) 

^  ^  8U  U\7rU  ^  mU^  1024t/‘’  '  ^ 

4  Curved- Surface  Specular  Region  Correction 

Details  of  the  curved-surface  specular  region  correction  are  presented  in  [22J.  Re¬ 
ferring  to  Figure  51,  the  term  of  the  diffraction  coefficient  which  “blows”  up  in  the 
specular  direction  is  corrected  by  the  inclusion  of  a  transition  function.  This  result 
goes  to  the  diffraction  point  solution  outside  the  caustic  region,  making  it  a  uniform 
solution.  The  specular  terms  from  the  two  diffraction  points  with  the  correction  are 


given  as 
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The  radius  of  curvature  of  the  diffractiug  edge  is  a, ,  and 


X  —  K  cot  a,.  Jtan^  a,  (cos^  v  +  c,^  sin^  +  ij  ^  sin  0„ 

2cm  j  cos  , 


(204) 


where  the  0’s  are  measured  from  the  curved  surface,  as  shown  in  Figure  51  and  e, 
is  the  ellipticity  of  the  curved  surface.  The  transition  function  is  given  as 


T{z)  =  T,{z)  =  -zM(h2.b,-z) 


(205) 


or 

T{z)^T2{z)  =  ~zU{\,2.^,-z\  (206) 

where 

f/(l,2.5,-2)  =  -?il/(l,2.5,-z)  +  ^(-z)-U-%  (207) 

and  M  and  U  are  confluent  hypergeometric  functions.  The  small  and  the  large 
argument  of  the  transition  function  are  the  following: 
for  \z\  <<  1, 


for  |z|  >>  1, 


/2z  2z'-'  z*\ 


TiU) 


(  2z  2z^ 

(v^  y  T  t) 


1 _ l_ 

sfz  2z  42^ 


r,(c)~,-±-±. 

The  new  diffraction  coefficient  can  be  expressed  as: 


=T 


D  {Qr)  =  ~0'0D^{Q..)  -  ^'^DUQ, ), 


(208) 

(209) 

(210) 
(211) 

(212) 
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where 


2n\/2nk 


:  sin  /?„  I  \ 


(4>  -  4>') 


R^j,  cot  ^ 


7r  +  (4>  +  ^')\ 
2n  j 


2n 


+  cot 


J)  H  CO. 


i>') 


2n 

T{±jXa,.) 


+ 


(213) 


R,j,  =  Tl. 


(214) 
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SECTION  5 


Top- Hat 

1  Introduction 

The  study  of  scattering  from  the  top-hat  has  been  of  interest  for  a  long  time.  Among 
some  of  the  latest  work  is  [26].  Which  uses  physical  optics  to  get  an  approximate 
backscatter  result.  In  order  to  get  a  more  accurate  UTD  result,  the  2-D  dihedral 
discussed  in  Section  3  can  be  extended  to  get  3-D  results  for  a  wide  variety  of  prob¬ 
lems.  The  solution  of  the  triangular  dihedral  in  the  principal  plane  is  discussed  in 
[13]  is  the  starting  point.  A  few  modifications  are  required  to  the  UTD  2-D  analysis 
of  the  dihedral  in  order  to  obtain  the  solution  to  the  top-hat.  The  high  frequency 
solution  is  validated  against  a  moment  method  body  of  revolution  result  [27].  For  a 
90°  dihedral,  as  seen  in  Section  2,  the  number  of  important  terms  decreases  signifi¬ 
cantly  for  the  backscatter  case.  It  is  shown  that  a  complete  solution  can  be  obtained 
by  including  the  D,  DR  and  RD,  DD,  and  DRD  terms. 

For  the  top-hat,  a  UTD  result  is  obtained  from  the  important  terms  which  is 
valid  in  all  regions  for  backscatter  and  is  reasonably  simple  to  program.  Also,  in 
the  dihedral  region,  the  dominant  result  is  further  simplified  to  provide  approximate 
results.  The  parameters  for  the  top-hat  are  shown  in  Figure  52. 
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Figure  52:  Top- hat  geometry. 
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2  Double  Diffraction  Simplified  to  Retain  the 
Significant  Part 


In  order  to  obtain  a  reiult  which  is  valid  in  all  regions  and  is  simple  to  program,  it  is 
necessary  to  extract  the  dominant  part  of  the  double  diffraction  for  the  backscatter 
case  when  the  two  diffractions  are  occurring  from  the  same  point.  This  is  only 
possible  in  terms  such  as  DRD. 

The  double  diffraction  formulation  for  the  joint  wedge,  discussed  in  Section  6,  is 
used  to  get  the  dominant  part  of  the  the  double  diffraction  coefficient.  In  this  case, 
m  =  n-2  =  n,  a,,  =  o,  =  a,  ‘t,,  -  and  d  are  the  same  as  defined  in  Section 

6  for  the  joint  wedge.  In  order  to  get  the  dominant  contribution,  it  is  necessary 
to  choose  #  =  TT  —  (mr  —  ^),  where  the  angle,  (^,  is  measured  from  the  ‘o’  face,  as 
shown  in  Figure  9.  Also,  the  distance  parameter,  d,  is  the  distance  traced  by  the 
ray  between  the  first  und  the  second  diffraction  points.  For  this  case,  the  double 
diffraction  coefficients,  P,  and  Pi,,  Equations  (45)  and  (46),  can  be  simplified  to 


F.(Qi,Qi)  =  P.(^,n,d) 


oTTjn  sin  (2,,) 


(215) 


for  a  TM  plane  wave  (with  the  electric  field  parallel  to  the  wedge),  and  for  a  TE 
plane  wave  (with  the  electric  field  perpendicular  to  the  wedge). 


(<?„<?,)  =  Ph  =  col'  (1) 


(216) 


For  the  top-hat,  these  equations  are  u'led  to  get  the  dominant  contribution  from 
D\RDi,  D:iRD:i  and  D^RDi  terms. 
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3  Modifications  Required 

Basically,  rim  caustic  correction  and  spread  factors  are  required  to  modify  the  2-D 
dihedral  solution  to  obtain  the  solution  for  the  top-hat.  In  Section  3,  it  is  shown  that, 
for  the  backscatter  case,  the  most  dominant  terms  for  the  90°  dihedral  are  diffrac¬ 
tion,  diffracted-reflected  (DR),  reflected-diffracted  (RD),  and  diffracted-reflected- 
diffracted  (DRD).  The  dominant  contribution  from  these  terms  are  used  to  get  the 
result  for  the  top-hat. 


3.1  D,  RD,  DR  and  DRD  terms 


For  the  backscatter  case,  referring  to  Figure  52,  using  Equation  (25),  the  diffracted 
terms  from  the  disk  is  given  by 
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n  -e 

jri’t  —  jn  . 


2\/27rfc 


1  -  R.. 


sin  (0) 


T<.  (2fc6sin  (0)) 


2  sin  (^) 
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=  U',J. 
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2\/2rk 
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sin  (ff) 


T,:(-2kbsm(0)) 


2  sin  (0) 


0  <  0  <  tan-'  [^]  , 


f  <  ^  <  TT, 


where 


=  < 


E,  E  field  perpendicular  to  the  pattern  plane 
H,  H  field  perpendicular  to  the  pattern  plane. 

Using  Equation  (13),  the  diffracted  terms  from  the  cylindrical  surface  are 
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„  2e-  '"/'sin 

_  r/<  .  V  »  / 

'-'ryl  —  '■■'f.h 


3v'^fc 


cos 


(?)  -  > 


(217) 


(218) 


(219) 
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T,  (2fcasiii  (^j) 


CEI 

y  2  sin  (6) 


^jlkh  <  >is(0). 


^  ^  _1  ^ 

0  <  -  +  tan  - -  , 

z  0  —  a 


(220) 


Op-Jir/l 


ze  ^  '  sin 


3\/2nk 


2fcasin(0)) 

cos(-^j  -cos(-LpjJ 

/  ^  <‘'s(fl)g- J  sin(0) . 


2  sin  (^) 


(221) 


The  diffracted- reflected  terms,  obtained  by  simplifying  Equations  (79)  and  (80) 
and  modifying  for  curvature  and  the  caustic  correction  are  given  as 

Ifiiih  ^  ^  {2kasm  (0)) 

\/27rk 


_ (t) 

3  [cos  (^')  -cos(^^)] 


+  j2kh  sin  (^) 


0  <  ^  <  tan' 


I  sin  (0)  ’ 
b  ~  a 

nr  ’ 


(222) 


f/"'"  -f  f/""'  =  Ul, 


2y/2nk 


2ft,,/. - 7^  -  i4fc  (6  -  a)  coc  (0) 

I  cos (^)  J 


2  sin  (0)’ 

<^<i- 


(223) 


For  the  DRD  terms,  as  discussed  in  Section  2,  only  the  dominant  part  is  retained. 
The  dominant  part  of  DRD,  where  the  diffractions  are  from  the  disk,  is 


(<!•  -  0  -  ^^-,n  2,rf  -  2(6  -  a)) 
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\  (pi  +  (6  -  a))  (p-^  +  (6  -  a)) 


P2  e 


,j2k(h-n)  sin(^7) 


where 


Pi 


P  == 


sin(^)  —  1  ’ 

((6  -  g)  +  pi)a 
a  -I  2((6  -  a)  +  p))’ 
6((6-a)  +  pO 


(224) 

(225) 

(226) 
(227) 


6  +  (sin(0)  -  1]  [(6  -  a)  +  p' j’ 
and  P,j,  ($,n,  <f)  is  given  by  Equations  (215)  and  (216).  The  dominant  part  of  DRD 

terms,  where  the  diffractions  are  from  the  cylindrical  surface,  are 

Ijihnih  ^  ($  =  e,n  =  2,d  =  2h)Tr(2kasin(0)) 


^  2(a  +  hsin(fl)) sin(0) 

0  <  0  < - tan"' 
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b  —  a 


(228) 


and 


^  ,  ^  ^  ^  ^  2,rf  =  2/i)  7;  (-2fcasin  (^)) 


V  2  (o  +  h  sin(0))  sin(0) 


g>  f  g  '  ”^(0)  ^J^kn  s;n(ff) . 


d<e  <  -.  (229) 

4  Dominant  Result  in  the  Double  Bounce  Re¬ 
gion 


It  is  shown  in  Section  2  that  the  diffracted- reflected  type  terms,  RD  and  DR,  may 
be  combined  to  give  a  dominant  result  in  the  double  bounce  region  for  the  dihedral. 
For  the  top-hat,  a  result  in  the  double  bounce  region  can  be  obtained  by  modifying 
Equations  (79)  and  (80), to  take  into  account  the  curvature,  so  that 


rr/^O.-i  +  Dfl  W 
^  Hnfvhun  ril 


d,.„,  <  6  <  tan ' ' 


b  —  a 
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(230) 
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and 


U. 


I),  li+iin, 

(iowinant 


=  Ulj, 


2a 


(6  — a)cos(0);  tan' 


<0  <  --Orn,.  (231) 


A  sin(0) 

In  Equations  (230)  and  (231),  is  an  angle  that  is  chosen  such  that  the  specular 
from  the  disk  and  the  top  of  the  cylindrical  rim  are  not  significant.  This  is  necessary 
since  these  equations  do  not  consider  singly  diffracted  fields.  A  reasonable  choice  for 
Oral  is  approximately  10°. 


5  PO  on  the  Axis  of  the  Top-Hat 

Along  the  axis  of  the  top-hat,  the  cotangent  terms  of  the  diffraction  coefficients  that 
tend  to  “blow”  up  and  combine  to  give  the  dominant  result.  Not  surprisingly,  the 
dominant  result  obtained  is  exactly  as  one  would  predict  from  physical  optics  [20]. 

As  expected,  at  the  z-axis,  the  reflection  boundary  cotangent  terms  associated 
with  the  diffraction  from  the  disk  (Equations  (217)  and  (218)),  give  the  area  of  the 
disk.  The  reflection  boundary  cotangent  terms  associated  with  the  diffraction  from 
the  upper  face  of  the  cylinder  (Equations  (220)  and  (221)),  give  the  area  of  the  end 
cap.  The  incident  boundary  terms  from  D:tR,  RD.)  and  DuRDi  (Equations  (222), 
(228)  and  (229)),  give  the  negative  area  of  the  end  cap.  Thus  the  apparent  “infinites” 
combine  to  give  a  result  which  equals  the  area  of  the  disk  minus  the  area  of  the  end 
cap  of  the  cylinder  plus  the  area  of  the  end  cap  of  the  cylinder  with  the  phase  taking 
into  account  twice  the  electrical  distance  of  the  cylinder.  The  RCS  along  the  axis  of 
the  top-hat  is  appropriately  given  as 

-V"  I 

6  Results 

Figures  53  and  54  compare  the  results  obtained  for  the  66  and  polarizations  by 
the  UTD  computer  code  and  the  Body  of  Revolution  code  at  10  GHz  for  a  top  hat 
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—  + 


—  ik2h 
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Figure  53:  Backscatter  field  with  66  polarized  field.  The  solid  line  is  UTD  solution 
and  the  daahed  line  is  the  BOR  solution. 


0  (deg.) 


Figure  54:  Backscatter  field  with  4>4>  polarized  field.  The  solid  line  is  UTD  solution 
and  the  dashed  line  is  the  BOR  solution. 
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of  cyJiiiHer  of  dimensions  o=  2  in,  h=  4  in,  and  h=  4  in.  The  results  compare  well 
with  the  body  of  revolution  code. 


7  B  ruder  hedral 

A  bruderhedral  is  a  reflector  that  is  an  angular  sector  of  a  top-hat.  (See  Figure  55.) 
Fundamentally,  the  solution  to  the  bruderhedral  is  the  same  as  that  of  the  top-hat, 
discussed  above.  In  the  region,  0  <  ^  <  90®,  the  scattered  fields  are  determined 
from  Equations  (217),  (220),  (222),  (223),  (224)  and  (228)  with  T,.  (U)  replaced  by 
Tf„riiai{U)  which  is  defined  as: 

T^rOai  (U)  =  -  {U)  )  r""”"  cos  (233) 

where  <  <l>  <  <l>bounj  is  the  angle  made  by  the  sector  of  the  bruderhedral. 

8  Results 

Figures  56  and  57  compare  the  results  obtained  for  the  90  and  <f><l>  polarizations 
by  the  UTD  computer  code  and  the  measurements  at  10  GHz  for  a  bruderhedral  of 
dimensions  a=  107  cm,  b=  115.9  cm,  c=  30.5  cm  and  h=  8.9  cm.  The  measurements 
are  courtesy  of  D.  Blejer  of  MIT,  Lincoln  Laboratory.  The  lip  of  the  model  used  in 
the  actual  measurements  is  believed  to  have  edges  parallel  to  the  y-axis,  as  shown  in 
Figure  55.  The  UTD  results  have  edges  along  constant  radial  lines.  The  differences 
are  expected  to  be  minimal  since  a  is  much  larger  than  the  difference  between  h  and 
o.  The  results  compare  well  with  the  measurements. 
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Figure  56:  Backscatter  field  with  00  polarized  field.  The  solid  line  is  UTD  solution 
and  the  dashed  line  is  the  measurement. 


Figure  57:  Backscatter  field  with  <t><f>  polarized  field.  The  solid  line  is  UTD  solution 
and  the  dashed  line  is  the  measurement. 
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SECTION  6 


Bi-Cone 

1  Introduction 

The  methods  used  to  analyze  the  top-hat  in  the  previous  section  are  very  similar 
to  those  required  to  obtain  the  scattered  fields  in  the  far-zone  from  the  bi-cone.  In 
general,  all  the  terms  discussed  for  the  UTD  2-D  solution  of  the  dihedral  in  Section  3 
have  to  be  modified  to  take  into  account  curvature  and  specular  caustics.  However, 
by  putting  a  restriction  on  the  bi-cone  that  the  sum  of  the  cone  angles,  a,  i  and 
a,  2  in  Figure  58,  is  a  right  angle,  for  the  backscatter  case,  the  number  of  important 
terms  decreases  significantly  as  seen  in  Section  6.  For  this  case,  a  dominant  UTD 
result  is  obtained  which  is  valid  in  all  regions  and  is  simple  to  program.  Also,  in  the 
dihedral  region,  the  dominant  r-isult  is  further  simplified  by  using  only  the  double 
bounce  terms  to  provide  approximate  results.  Figure  58  shows  the  parameters  for 
the  bi-cone. 

2  Modifications  Required 

Basically,  it  is  only  necessary  to  change  the  caustic  corrections  and  spread  factors 
in  the  2-D  dihedral  solution  to  obtain  a  result  for  the  bi-cone.  Two  types  of  caustic 
corrections  are  needed.  One  is  the  rim  caustic  which  is  discussed  in  Section  3  and 
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Figure  58:  Bi-Cone  geometry. 
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used  in  the  solution  of  the  top-hat.  in  Section  5,  anrl  the  other  is  the  curved-surface 
specular  region  correction  whicl,  is  examined  in  Section  4. 


3  Dominant  Result 


The  bi-cont  has  the  same  type  of  dominant  mechanisms  as  the  2-D  dihedral.  In  order 
to  simplify  the  problem,  it  is  necessary  to  choose  cone  angles  such  that  a,  i  -1-  a,  2  = 
90°.  Thus,  as  for  the  2-D  dihedral  discussed  in  Section  6,  the  dominant  result  is 
obtained  from  diffracted  (D),  RD,  RD  and  DRD  terms  for  the  bi-cone. 


3.1  D,  RD,  DR  and  DRD  terms 


The  expressions  shown  below  are  for  backscatter  from  a  bi-cone  of  parameters  given 
in  Figure  58.  Also  note  that 
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Using  Equations  (25)  and  (26),  the  diffracted  fields  from  the  upper  rim  of  the 
bi-cone  of  radius  c  is  given  by 
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The  diffra  cted  tiei  J  from  the  lower  rim  of  the  bi-cone  of  radius  b  is  given  by 


(240) 


(241) 


The  purpose  of  the  curved  surface  specular  caustic  correction  for  the  DR  and 
DRD  terms  is  twofold:  to  ensure  continuity  of  the  fields,  and  to  obtain  the  physical 
optics  result  for  the  visible  part  of  the  cone  in  the  specular  direction.  This  requires 
that  for  the  diffraction-reflection  type  terms,  the  caustic  correction  is  associated 
with  the  incident  shadow  boundary.  After  simplifying  Equations  (70)  and  (71)  and 
modifying  for  curvature  and  caustic,  the  diffraction-reflection  terms  are  given  as 
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and 
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In  the  above  equation  for  Di  R  and  RDt  terms 
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and 
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For  the  DRD  terms,  as  discussed  in  Section  2,  only  the  dominant  part  is  retained 
for  the  joined  wedges  solution.  The  dominant  part  of  the  DRD  term  in  the  region 
where  z  is  positive  (especially  around  6  =  0,2)  is 
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In  Equation  (252) 
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4  Dominant  Result  in  the  Double  Bounce  Re¬ 


gion 


It  is  shown  in  Section  2  that  the  diffracted-reflected  type  terms,  RD  and  DR,  may 
be  combined  to  give  a  dominant  result  in  the  double  bounce  region  for  the  dihedrad. 
Modifying  Equations  (79)  and  (80),  to  take  into  account  the  curvature,  a  result  in 
the  dihedral  region  for  the  bi-cone  can  be  given  as 


fjnn,+n,ii  ^  fj, 


sin  (0  -  0,  2) ; 


\/A  \  {pV  +  ^i) 

a,  2  +  Q  vut  <e  <  a<,2  +  tan  ' 


(260) 


and 


rjlhH+Hlh  _  jj,  , 

^dominant  ^ /Y  \ 


q,2  +  tan 


-I 


sin  {6  +  Ori) ; 
<6  <Tf  -  arA  -  Or,,,. 


ip?  +  oi) 

lo 


(261) 


Using  Equations  (244)  to  (247)  in  Equation  (260),  and  using  Equations  (248)  to  (251) 
in  Equation  (261),  the  dominant  contribution  from  the  double  bounce  is  obtained. 
An  angle,  6r„i,  is  chosen  so  that  the  specular  from  the  curved  sides  of  the  cones 
are  not  significant.  This  is  necessary  since  these  equations  do  not  consider  singly 
diffracted  fields.  A  reasonable  choice  for  Oru,  is  approximately  10°. 


5  Results 

Figures  59  and  60  show  the  results  obtained  for  the  90  and  <{>(p  polarizations  by  the 
UTD  computer  code  and  the  moment  method  body  of  revolution  code  [27]  at  10  GHz 
for  bi-cone  of  a  =  0.5  in  (0.42  A),  6  =  c  =  2.0  in  (1.69  A),  and  l\  =  I2  ~  1.5\/2  in 
(1.8  A).  Considering  the  small  size  of  the  bi-cone,  the  results  compare  well  with  the 
body  of  revolution  moment  method  code. 


100 


'  0  20  40  60  80  100  120  HO  160  60 

0  (deg.) 


Figure  59;  Backscatter  field  with  66  polarized  field.  The  solid  line  is  UTD  solution 
and  the  dashed  line  is  the  BOR  solution. 


Figure  60:  Backscatter  field  with  4><i>  polarized  field. The  solid  line  is  UTD  solution 
and  the  dashed  line  is  the  BOR  solution. 


SECTION  7 


Parallel  Plate  Cavity  and  Circular 
Cavity 

1  Introduction 

The  study  of  scattering  from  a  parallel  plate  cavity  and  its  three-dimensional  coun¬ 
terpart,  a  circular  cavity,  has  been  of  much  interest  lately  [28,  29).  These  shapes 
may  be  used  to  model  many  common  objects.  For  example,  the  parallel  plate  cavity 
can  be  used  to  represent  the  dominant  scattering  from  a  truck  bed  and  the  shallow 
circular  cylinder  cavity  resembling  a  cake  pan  can  be  used  to  represent  the  dominant 
scattering  from  a  wheel  hub,  cannon  barrels,  jet  inlet  structure,  et  cetera.  The  the¬ 
ory  and  the  results  discussed  in  this  section  are  valid  for  cavities  of  arbitrary  depths. 
However,  since  the  terms  in  the  non  principal  plane  are  not  included,  the  smaller 
the  height  to  width  ratio  of  the  cavity,  the  wider  the  range  of  angles  over  which 
the  result  is  valid.  Also,  the  dominant  terms  of  a  shallow  cavity  are  the  same  over 
much  of  the  region  (the  smaller  the  height  to  width  ratio  of  the  cavity,  the  larger  is 
the  region)  as  for  a  dihedral,  which  has  been  studied  in  detail  in  Section  3.  In  this 
section,  three  dots,  “•  ■  will  represent  one  or  more  reflections  within  the  cavity. 

2  Dominant  Terms  for  the  Parallel  Plate  Cavity 

Consider  the  two-dimensional  geometry  shown  in  Figure  61.  For  the  backscatter 


z 

^  B 


Figure  61:  Two-dimensional  parallel  plate  cavity  geometry. 
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case,  looking  along  the  z-axis  {$  =  0"),  the  important  terms  are  the  difTrartions 
(Di  and  D2),  diffracted-reflected  {DiRf,  and  Ri,Di)  and  the  diffracted-  reflected- 
diffracted  {D]Ri,Di  and  D^RiiD^)  terms.  See  Figure  62.  Ri,  refers  to  the  reflection 
from  the  base  of  the  cavity.  From  the  study  of  the  basic  dihedral  mechanisms  in 
Section  2,  it  is  expected  that  as  d  departs  from  0°,  the  diffracted-reflected  {D\Ri, 
and  Ri,D\)  terms  will  dominate  the  solution.  However,  the  angular  range  over  which 
the  diffracted-reflected  terms  dominate,  will  depend  on  the  ratio  of  the  height  to  the 
width  of  the  cavity.  The  shallower  and  wider  the  cavity  the  larger  is  the  angular 


range. 

As  0  increases,  one  set  of  terms  get  shadowed  and  another  set  of  higlier  order 
terms  appear.  For  example,  once  the  DR  and  RD  terms  get  shadowed  the  DRR  and 
RRD  terms  appear,  upon  shadowing  of  DRR  and  RRD  terms,  DRRR  and  RRRD 
terms  appear,  et  cetera.  From  a  geometrical  point  of  view,  the  D)  R  •  •  •  and  •  •  •  RDi 
terms  are  visible  in  the  region  (•  •  •  represent  reflections) 


<  0  <  tan  ' 


(2m  -f-  1)  d 
^  ’ 


(262) 


where  m  is  the  number  of  reflections  before  reflecting  from  the  base  of  cavity. 
For  example,  m=0  represents  D\Rt„  m=l  represents  DiRR4,R,  m=2  represents 
D\  RRRi,RR  and  so  on.  Rf,  refers  to  the  reflection  from  the  base.  (See  Figure  62.) 
Also,  the  D'xR  •  •  •  and  •  •  •  RD2  terms  are  visible  in  the  region 


where  m  is  the  number  of  reflections  before  reflecting  from  the  base  of  cavity. 
For  example,  m=0  represents  D|R(,R,  m=l  represents  D\  RRi,RR,  Tn—2  represents 
D\RRR\,RRR  and  so  on.  Again,  R^  refers  to  the  reflection  from  the  base  (see  Fig¬ 
ure  62).  Given  the  angle  6,  the  value  of  m  (number  of  reflections  before  reflecting 
from  the  base  of  the  cavity)  may  be  determined  using  Equations  (262)  and  (263). 
Rearranging  Equation  (262), 

2hls.n6  —  d  hlanO  , 

-IT-  <  ”■  <  — 
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and  if  tb<*  range  of  m  inrlndes  an  integer,  then  for  the  given  angle,  0,  the  /)(•••  term 
is  present.  Thus,  the  diffraction  occurs  from  the  upper  left  edge  of  the  parallel  plate 
cavity,  and  the  integral  value  of  m  is  the  number  of  reflections  before  reflecting  from 
the  base  of  the  cavity.  Rearranging  Equation  (263), 

h  tan  S  —  d  2h  tan  9  —  d  ,  . 

—1—  ^  ^  -  hi  - 

and  if  the  range  of  m  includes  an  integer,  then  for  the  given  angle,  0,  the  D2  •  •  •  term 
is  present.  Thus,  the  diffraction  occurs  from  the  upper  right  edge  of  the  parallel 
plate  cavity,  and  the  integral  value  of  m  is  the  number  of  reflections  before  reflecting 
from  the  base  of  the  cavity.  Note,  in  general,  for  a  particular  6  either  Equation 
(264)  or  (265)  will  have  an  integral  value  in  its  range.  If  both  the  equations  have  an 
integral  value  in  the  range,  then  there  is  a  shadow  boundary  at  that  angle. 


2.1  Di**-and  •••JDi  Terms 

Consider  the  case  of  diffraction  from  the  top  left  edge  (D| )  shown  in  Figure  63.  The 
field  from  Di  •  •  •  for  a  bistatic  angle  of  c,  where  e  is  a  small  number,  valid  in  the 
region  defined  by  Equation  (262),  is  given  as 


U 


n(2"'  +  2)  . 


-  (*  +  »')  j 


,>*  (n  r'+rj  -r) 


(266) 


where  m  is  the  number  of  reflections  before  reflecting  from  the  base  of  the  parallel 
plate  cavity,  and 

f  =  sin  9  X  +  cos  9z,  (267) 

r  =  sin  (^  —  e)  z  +  cos  {9  —  e)  z,  (268) 

r,  =  —a  X  +  h  z  (269) 


is  the  position  vector  of  the  diffraction  point  (/?)), 


=  —  (4m  +  2)a  X  —  h  z 


(270) 
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is  the  position  vector  of  the  diffraction  point  {D\)  imaged  through  all  reflection, 


and 


(271) 


(272) 


^'  =  n-e. 

Using  Equations  (267)  to  (272)  in  Equation  (266), 

JjHf-  ^  ^  ^  lm+2)n  siiiflgjtf  |(  lrii+ l)a  <:i>sO-/i  sinO]  f273) 

v®  ■  K  ) 

The  field  from  ■  •  •  Hi,  valid  in  the  region  defined  by  Equation  (262),  is  given  as 


U- 


■1), 


iiir”  cot '''-<*  +  * 


4i\/2irk 


(274) 


with  Equations  (267)  to  (270)  the  same  as  above  and 


$  =  rr  —  (0  —  e) 


and 


(275) 


(276) 


Equation  (274)  may  now  be  simplified  to 

i,  t'. 

\e/ 

Thus,  the  total  field  from  the  diffraction  occurring  from  the  terms  with  Z)*,  for  the 


g  lni+2)o  sin  0^jkf[n  mnO-kh  sin  B] 


(277) 


backscatter  case,  (Equations  (273)+(277)  in  the  limit  as  e  — »  0),  valid  in  the  region 
defined  by  Equation  (262),  is 


IJtnlal  Irjl  ^  g->*(  lm+2)' 


2.2  Po  •  •  •  and  •  •  •  Z??  Terms 


[h  sin  9  —  2ma  cos  9] . 


(278) 


Consider  the  terms  with  diffraction  from  the  top  right  edge  (H^)  shown  in  Figure 
64.  The  field  from  Z)2  •  •  •  for  a  bistatic  angle  of  e,  where  e  is  a  small  number,  valid 
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in  the  region  defined  by  Equation  (263),  is  given  as 


I/"’  ~  _ 


e 

4\/2irk 


cot  ^  j  (279) 


where  m  is  the  number  of  reflections  before  reflecting  from  the  base  of  the  parallel 
plate  cavity  and  Equations  (267)  and  (268)  are  the  same  as  before,  and 


Fj  =  a  X  -i-  h  z 


(280) 


is  the  position  vector  of  the  diffraction  point  (D-z), 


f"  =  —  (4m  +  3)  a  i  —  A  z 


(281) 


is  the  position  vector  of  the  diffraction  point  (D^)  imaged  through  all  reflection. 


$  =  (^-e). 


(282) 


and 

$'  =  IT  +  e.  (283) 

Equation  (279)  can  now  be  simplified  to 

(jfh  ■  ^  ^  g-iMlm+2)rt*iiiflgj>,[(lm+:«)oc.H>0-ft*infl]  (284) 

v^Vc/ 

The  field  from  •  •  •  D^,  valid  in  the  region  defined  by  Equation  (263),  is  given  as 

U  ih  ^  cot  ( ~  .  (285) 

4V27rfc  ’  \  4  / 

Using  Equations  (267),  (268),  (280),  (281), 


#  =  jr  +  (^  -  e). 


and 


=  9. 


Equation  (285)  may  be  simplified  to 


e  /I 


y/2vk 


^  —  lm  +  2)n  ^in  I'*  <  i*'"  ^1 


(286) 


(287) 


(288) 


110 


The  total  field  from  the  diffraction  occurring  from  the  terms  with  /)-,  for  the 
backscatter  case,  (Equations  (284)+(288)  in  the  limit  as  e  — >  0),  valid  in  the  re¬ 
gion  defined  by  Equation  (263),  is 


i/? 


U,oiair,„i.i  _  (2(m  l)acos0  -  hsme] .  (289) 


It  can  easily  be  shown  that  when  both  the  diffractions  of  Di  •  •  •  and  D2  •  ■  •  are  at 
the  shadow  boundaries  defined  by  Equations  (262)  and  (263),  the  terms  Di  •  •  •  and 
D2  •  •  •  are  continuous,  and  Equations  (278)  and  (289)  may  be  further  simplified.  At 
the  boundary. 


0  =  61  =  cot  * 


2h 

(2m  +  l)d 


m  =  0,1,--  -  , 


(290) 


the  field  from  D  ■  •  •  R  type  terms  is 


=  -7=a  cos 

V  A 


(291) 


At  the  shadow  boundary, 

0  =  0,!  —  cot~' 
the  field  from  D  •  •  ■  R  type  terms  is 


(m  f  1)  d 


;  m  =  0,1,---, 


(292) 


U 


n  +  n 


lo^o,.  =  0. 


(293) 


From  Equations  (264)  and  (265),  as  0  approaches  90°,  more  and  more  terms  are 
required.  However,  from  Equation  (291)  it  is  observed  that  as  0  increases  towards 
90°  the  result  approaches  zero,  so  it  is  no  longer  necessary  to  include  these  higher 
order  terms.  This  should  be  expected  physically  since  less  and  less  of  the  cavity  is 
visible  and  the  broadside  reflection  from  the  outer  cylinder  becomes  dominant. 


2.3  D '  •  •  D  Terms 

Also  of  importance  are  the  D,  --Di  and  Dj  -  D^  terms.  Although,  for  Ihe 
backscatter  case,  these  terms  do  not  contribute  to  the  continuity  of  the  fields,  they 
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are  as  import.anf  as  <lie  f)  ■  ■  ■  and  P  type  terms.  Tt  is  fonnH  that  ilie  dominant 
part  of  the  terms,  D\  •  •  •  D-i  and  D^--  •  D\^  associated  with  both  the  diffractions  be¬ 
ing  near  the  shadow  boundaries,  cancels.  However,  a  significant  part  remains  from 
the  non-dominant  part  of  the  terms,  D\  ■  •  ■  D-2  and  D2  ■  •  ■  D],  associated  with  only 
one  of  the  diffractions  being  on  the  shadow  boundary. 

Typical  results  from  terms  of  the  type  D  ■  ■  ■  and  •  •  ■  D  (Equations  (278)  and  (289) 
valid  in  the  regions  given  by  Equations  (262)  and  (263),  respectively)  are  shown  in 
Figure  65.  At  the  boundaries.  Equations  (290)  and  (292),  the  D  ■  and  •  ■  •  D  terms 
go  through  peaks  (Equation  (291))  and  nulls  (Equation  (293)),  respectively.  In 
general,  as  shown  in  Figure  65,  the  D  ■  •  •  D  terms  relatively  smooth  out  the  peaks 
and  the  nulls  of  the  D  •  •  ■  and  •  •  ■  D  terms. 

3  Parallel  Plate  Cavity  Verses  Dihedral 

From  the  dihedral  study  in  Section  2,  and  Figures  15  and  16,  it  is  observed  that  the 
DD  across  the  face  of  the  dihedral  balances  the  slope  discontinuity  of  the  DR  and 
RD  terms  for  the  backseat  ter  case.  However,  as  seen  in  the  above  section,  for  the 
parallel  plate  cavity,  the  D  •  •  •  D  terms  relatively  smooth  out  the  ripple  due  to  the 
D  -  •  •  and  •  •  D  terms.  Thus,  the  response  of  a  parallel  plate  cavity  is  not  as  sharp 
as  compared  to  that  of  the  dihedral. 

4  Results 

Figures  66  and  67  compare  the  results  obtained  by  including  terms  up  to  12  reflec¬ 
tions  after  diffraction  and  the  Method  of  Moments  at  10  GHz.  The  dimensions  of 
the  cavity  are  d=4  in  (3.387  A)  and  h=8  in  (6.773  A).  As  expected,  the  results 
compare  well  with  the  method  of  moments. 
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Figure  65:  Importance  of  D---D  terms  for  a  parallel  cavity  d=4  in  and  h=8  in 
with  <()<f>  polarized  field. 
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5  Circular  Cavity 


Consider  the  circular  cavity  shown  in  Figure  68.  Apart  from  the  terms  present  in 
the  parallel  plate  cavity,  terms  with  single  and  double  diffractions  and  reflections 
are  present  that  are  not  in  the  principal  plane  of  the  cavity.  However,  terms  with 
diffraction  points  not  in  the  principal  plane  are  not  considered  in  this  study.  It  is 
found  that  for  circular  cavities,  the  smaller  the  height  to  width  ratio,  the  larger  is 
the  region  over  which  the  terms  in  the  non-principal  plane  are  insigniflcant.  The 
circular  cavity  result  is  obtained  by  modifying  the  terms  of  the  parallel  plate  cavity 
by  the  spread  factor  and  the  caustic  correction.  The  caustic  corrections  are  basically 
the  same  as  that  for  the  top-hat  discussed  in  detail  in  Section  5. 


6  Spread  Factor  Modification 

In  order  to  account  for  the  concave  geometry,  the  spread  factor  has  to  be  modified 
to  take  into  consideration  each  bounce  on  the  circular  cavity. 


6.1  Spread  Factor  for  R'  • '  D2 


Referring  to  Figure  69,  the  spread  factor  for  R  -  •  -  D2  will  be  derived.  For  the 
backscatter  case,  the  spread  factor  for  R  -  ■  ■  Df  is  the  same  as  for  Di  •  •  •.  In  general, 
for  m  reflections  after  reflecting  from  the  base  of  the  cavity,  the  spread  factor  may 
be  represented  as 


where 


P^prrnH  farinr 


_rl  _rl 

Pi  P? 


„t2 

P\  P2 


[P\'‘  +  5)  (P?'  +  J*)  {P?  +  5)  [P?  +  3) 

^r(2?ti+  I  I ) 

"  {p^'Us)  + 


P\  =  P\  - 


P2' 


r(2m+l) 

^  Pi 


=  00, 


-a 

2sin 


(294) 


(295) 

(296) 
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z 


•  •  /?2  type  terms  in  circular  cavity. 
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and 


s  = 


2a 

sin  0' 


(297) 


P2^ 


S  +  P2 


is  the  distance  between  two  points  of  contact  on  the  circular  part  of  the  cavity.  For 
the  second  reflection, 

(298) 

(299) 


Pi  — 


+7>5’)’ 

1 


1_2^(,  +  ^H) 


and  for  the  n  th  reflection 


P2 


P2  = 


1 

1  -  ^  (> + ft'""’) ' 


and  for  the  diffraction, 


r(2m+l)  , 

p  =  Pi'  '  +  S. 

Using  Equations  (295)  to  (302)  in  Equation  (294),  the  spread  factor  becomes 


r2«  r3-  r(2m+l)-l? 


_  f  rl  .  - 

ppprraH  factor  ■“  \p2  P2  Pi  *  *Pl 

6.2  Spread  Factor  for  R  -  ••  Di 


r- 


(300) 

(301) 

(302) 

(303) 


Referring  to  Figure  70,  the  spread  factor  for  R’--D\  will  be  derived.  For  the 
backscatter  case,  the  spread  factor  for  iZ  •  •  -  Z?2  is  the  same  as  for  D2  •  •  ••  In  general, 
for  m  reflections  after  reflecting  from  the  base  of  the  cavity,  the  spread  factor  may 
be  represented  as 


P spread  factor 


A(pr,'+,)(p',^/j)(pr,^  +  3)  ' 

_  1 

■■(p^^"’  +  s)(p5'^"'  +  s)^J  ’ 


(304) 


where 


Pi  =  Pi 


=  Pi 


=  00, 


(305) 
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•  D\  type  terms  in  circular  cavity. 
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rl 


Pi  = 


and 


—  a 

2  sin  9’ 
2a 


s  = 


sin  9’ 


(306) 


(307) 


is  the  distance  between  two  contact  points  in  the  circular  part  of  the  cavity.  For  the 
second  reflection, 


Pi""  = 


Pi 


S  +  Pl^ 


i-^i^+Pi'y 

1 


and  for  the  n  tk  reflection, 


Pi  = 


Pi 


»  +  p2 


and  for  the  diffraction, 


and 


P  = 


1 

jp?”'  +  s) 

1 


(P2^"'  +  S) 


* 


Using  Equations  (305)  to  (313)  in  Equation  (304),  the  spread  factor  becomes 


P.^prf:ad  fnvtor  —  ^2  P2  P2  '  *P2  P 


However,  if  m  =  0  (i.e.,  for  the  terms  D\Rt,  and  /?(,£>]) 

[i 

—a  1 

2^9.  • 


(308) 

(309) 

(310) 

(311) 

(312) 

(313) 

(314) 


(315) 


6.3  Spread  Factor  for  D\”  'D2 


Referring  to  Figure  71,  the  spread  factor  for  Dy-’-Dj  will  be  derived.  For  the 
backscatter  case,  the  spread  factor  for  D\  •••Di  is  the  same  as  for  Di-'-Dy.  In 
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general,  for  m  reflections  after  reflecting  from  the  base  of  the  cavity,  the  spread 
factor  may  be  represented  as 


P^jyrrnH  fai-tiir  — 


Pi 


P\'P2 


P?P? 


» (pi  +  j)  (p’l '  +  s)  ip2  +  3)  (p’,^  +  a)  (p2^  +  a) 
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,  i 
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(316) 


where 


Pi  = 

—a 

(317) 

(sin^  +  sin)9)’ 

2a 

sin/3’ 

(318) 

P?  =  3, 

(319) 

Pi*  =  2s, 

(320) 

p\^”'  =  2ms, 

(321) 

for  the  first  reflection. 

.rj  _ 

P2  — 

pi  +  s 

(322) 

l-^(p.+o)’ 

-rl*  _ 

Pi  - 

1 

(323) 

for  the  second  reflection. 

p7  = 

p5‘  +  ^ 

(324) 

l-^W+o)’ 

pT  = 

1 

(325) 

and  for  the  n  th  reflection. 
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(326) 
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P2 
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and 


P2  — 


1  _  (pr2m  ^  j)' 

Using  Equations  (317)  to  (329)  in  Equation  (316),  the  spread  factor  is 


(329) 


Pikprrad  jnrior  — 


_  .rl  _r2»*r3» 

P^P2  P2  P2 


P?'"' 


(2m  +  1) 


(330) 


7  Cake  Pan  Results 


In  this  section,  the  results  obtained  by  the  UTD  method  mentioned  above  are  com¬ 
pared  with  the  modal  theory.  The  results  for  the  modal  theory  are  obtained  from 
[29].  Figures  72  and  73  are  for  backscatter  with  4><i>  and  09  polarizations  for  a  cake 
pan  of  both  the  radius  and  the  height  of  two  wavelengths.  Figures  74  and  75  are  for 
backscatter  with  and  00  polarizations  for  a  cake  pan  of  radius  two  wavelengths 
and  a  height  of  four  wavelengths.  Figures  76  and  77  are  for  backscatter  with  <i>i> 
and  00  polarizations  for  a  cake  pan  of  radius  two  wavelengths  and  a  height  of  ten 
wavelengths.  Figures  78  and  79  are  for  backscatter  with  and  00  polarizations  for 
a  cake  pan  of  radius  five  wavelengths  and  a  height  of  ten  wavelengths. 

From  a  simple  ray  picture  it  may  be  shown  that  the  first  D  ••  •  term  which  is  not 
in  the  principal  plane  and  hence  not  considered  in  this  solution,  appears  at 

0  =  tan~' 

Note  that  this  is  a  boundary,  and  the  solution  may  begin  to  deteriorate  even  before 
it.  This  corresponds  to  0  ~  54°  for  Figures  80  and  81,  ^  ~  35°  for  Figures  72  and 
73,  0  ~  19°  for  Figures  74,  75,  78  and  79,  and  ^  ~  8°  for  Figure  76  and  77.  Taking 
this  into  consideration,  it  is  observed  that  the  results  compare  well  with  the  modal 
theory. 


hy/2 


(331) 


oo 

m 

o 


c> 

<c 


Figure  72:  Backscatter  field  for  a  circular  cylinder  with  radius  2A  and  height  2A  with 
polarized  field  (UTD  solution  valid  for  0  <  35°). 


Figure  73:  Backscatter  field  for  a  circular  cylinder  with  radius  2A  and  height  2A  with 
polarized  field  (UTD  solution  valid  for  <  35°). 


124 


ANGLE  (DEG.) 


Figure  74;  Backscatter  field  for  a  circular  cylinder  with  radius  2A  and  height  4A  with 
<f>(j>  polarized  field  (UTD  solution  valid  for  B  <  19°). 


O  5  10  19  20  25  50  35  *0  45  SO 


ANGLE  (DEG.) 

Figure  75;  Backscatter  field  for  a  circular  cylinder  with  radius  2A  and  height  4A  with 
6B  polarized  field  (UTD  solution  valid  for  9  <  19°). 
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Figure  76;  Backscatter  field  for  a  circular  cylinder  with  radius  2A  and  height  lOA 
with  <l>4>  polarized  field  (UTD  solution  valid  for  9  <  8°). 
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Figure  77:  Backscatter  field  for  a  circular  cylinder  with  radius  2A  and  height  lOA 
with  99  polarized  field  (UTD  solution  valid  for  9  <  8°). 
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Figure  78;  Backscatter  field  for  a  circular  cylinder  with  radius  5A  and  height  lOA 
with  <f><f>  polarized  field  (UTD  solution  valid  for  0  <  19°). 


Figure  79;  Backscatter  field  for  a  circular  cylinder  with  radius  5A  and  height  lOA 
with  90  polarized  field  (UTD  solution  valid  for  9  <  19°). 
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Figure  80:  Backscatter  field  for  a  circular  cylinder  with  radius  4A  and  height  2X  with 
<i>4>  polarized  field  (UTD  solution  valid  for  6  <  54°). 


Figure  81:  Backscatter  field  for  a  circular  cylinder  with  radius  4A  and  height  2X  with 
Off  polarized  field  (UTD  solution  valid  for  0  <  54°). 
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SECTION  8 


Donut 

1  Introduction 

First  and  second  order  Geometrical  Optics  solution  and  their  caustic  corrections  for 
the  donut  are  obtained  in  this  section.  This  is  a  special  case  of  scattering  from 
elliptical  cylinders  discussed  in  detail  in  (30].  This  configuration  is  of  interest  since 
it  can  represent  tires  on  a  vehicle  and  with  Section  7  a  wheel  hub.  It  is  observed  in 
[30]  that  a  fourth  order  equation  has  to  be  solved  in  order  to  determine  the  reflection 
points  from  such  a  geometry.  However,  working  the  other  way  around,  with  one  of 
the  reflection  points  assumed  known,  a  simple  solvable  equation  can  be  obtained  to 
get  the  angle  of  incidence.  Interpolation  can  then  be  used  to  determine  the  reflection 
points  for  a  given  angle  of  incidence.  In  Section  2,  for  a  general  bistatic  case,  the 
points  of  reflections  for  the  double  bounce  are  determined. 

The  rest  of  the  section  deals  with  backscatter.  First,  the  GO  term  for  the  2-d 
cylinder-cylinder  doubly  reflected  term  is  obtained  in  Section  3.  This  is  used  later 
to  get  GO  field  for  the  donut.  The  first  order  GO  solution  for  the  donut  and  its 
caustic  correction  are  obtained  in  Sections  4  and  5,  respectively.  In  Section  6,  the 
second  order  GO  solution  for  the  donut  and  its  caustic  correction  is  discussed. 
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2  Cylinder- Cylinder  Interaction 

Consider  the  two-dimensional  geometry  shown  in  the  Figure  82.  In  order  to  compute 
the  spread  factors,  it  is  necessary  to  determine  the  points  of  reflection.  It  was  found 
in  Section  3,  that  for  the  double  bounce  to  occur,  the  angle  of  the  dihedral  must 
be  a  =  I  -I-  or  a  =  I  —  |^^|.  Thus,  for  a  given  bistatic  angle,  two  dihedral 

angles  are  possible  which  will  result  in  a  double  bounce.  Of  the  two  double  bounces, 
one  is  for  the  case  when  the  first  reflection  occurs  from  cylinder  of  radius  Oi,  R1R21 
and  the  other  is  for  the  case  when  the  first  reflection  occurs  from  cylinder  of  radius 
02,  R^Rx  (/2|  and  i?2  refers  to  reflections  from  the  cylinders  of  radius  and  02, 
respectively).  However,  for  the  backscatter  case,  there  is  only  one  dihedral  angle, 
a  =  |,  and  the  two  ray  paths,  RiRj  and  R^Rx-,  are  the  same.  We  will  determine  the 


Figure  82:  Cylinder  Cylinder  interaction. 
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reflection  points  for  one  of  the  rays,  R\R-i  or  R/Rt,  and  use  that  result  to  determine 
the  reflection  points  for  the  other  ray  path.  First,  the  ray  path  that  results  is  the 
dihedral  angle  being  greater  then  ^  will  be  considered;  i.e.,  a  =  |  +  |^^|-  This 
choice  is  necessary  in  order  to  avoid  the  ambiguity  in  the  relationship  between  sides 
and  angles,  the  sine  rule. 

From  Figure  82,  it  is  observed  that  if  0  >  6\  a  is  greater  then  ~  (i-e.,  a  = 
2  determine  the  angle  ^  for  R1R2.  However,  if  ^  <  0',  we  will 

determine  the  angle  0  for  R2R{.  From  Figure  82,  one  gets 


7 


(332) 


Using  the  sine  rule  on  the  triangle  formed  by  ni. 


h.2  and  the  y-axis. 


dt  d^ 


S\ 


S2 


sin(7r  —  a)  sin(|  —  7)  sin/3’ 


(333) 


where  s\  and  ^2  are  the  length  made  by  sides  ni  and  1x2,  respectively.  Using  Equation 
(333)  and  applying  sine  rule  on  the  triangle  formed  by  nj,  n2  and  the  two  points  of 
reflection. 


n  = 


(di  +  d2)cos7  —  ai  sin  a  8in(^'  +  7) 


(di  +  ^2)  sin  /3  —  02  sin  a  cos  {0  +  0) 
Rearranging  Equation  (334)  and  solving  for  the  incident  angle,  0, 

n  cos  0  —  cos{0  —  a  —  ^b)l 


0  =  tan 


-I 


where  0(,  is  the  bistatic  angle, 


[  n  sin;3  —  sin(/3  —  a  —  ^b)  J 


ft,  =  «  -  0'. 


(334) 


(335) 


(336) 


Equation  (335)  may  be  used  to  obtain  the  incident  angle,  0,  if  the  point  of  first 
reflection,  or  0,  is  known.  However  if  0  is  known,  which  is  normally  the  case,  then 
0  may  be  found  by  trial  and  error.  A  more  reasonable  approach  may  be  to  find  ^’s 
for  particular  values  of  0,  using  Equation  (335),  and  then  interpolating  to  find  0  for 
a  particular  0.  Now,  7  may  be  found  from  Equation  (332). 
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Now,  ronsider  the  other  rase,  when  nt  =  ^  —  j^^ji  referring  to  the  other  »IonI)le. 
bounce,  it  can  be  shown  that 


~ 


(337) 


and 


1/olhcr  — 


0-0' 

2 


(338) 


Thus,  one  double  bounce  occurs  at  reflection  points  defined  by  angles  jS  and  7,  and 
the  other  double  bounce  occurs  at  the  reflection  points  defined  by  angles  r  and 


'Jothcr  • 


3  Two-Dimensional  Cylinder- Cylinder  Doubly 
Reflected  Term 


Referring  to  the  two-dimensional  geometry  shown  in  Figure  83.  The  backscatter 
doubly  reflected  term  R\R2,  using  Equation  (1),  is  given  as 


(339) 


Ri  and  R2  are  defined  in  Equation  (1)  at  the  first  and  second  reflection  points, 
respectively,  and  p\'  and  p]’ ,  the  caustic  distances  for  the  first  and  the  second 
reflected  rays,  are  given  as 

_  osin(fl  -f  0) 


and 


J_  -  1  2 

s  +  pV  ^  a  cos  (0-1- 7) 


(341) 
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Figure  83:  Two-dimensional  Cylinder-Cylinder  showing  R1R2  interaction. 
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4  Donut;  First  Order  Solution 


Consider  modifying  the  2-D  double  cylinder  to  a  donut,  shown  in  Figure  84.  For  the 
backscatter  case,  the  two  single  reflected  fields,  from  Equation  (1),  are  given  as 


£'  =  t  (0)  • 


(342) 


where  p\  and  obtained  using  Equation  (5).  The  radii  of  curvatures  are  R\  = 

T^slud  ®  ^2  =  <!•  Also,  using  ^’  =  0,  ^  and  ^2  =  |  in  Equation  (5),  the 

caustic  distances  become 

(343) 


Pi  o 


=  T 


2’ 

d  a 
4sin^  ^  2 


(344) 


In  Equation  (342),  fv  =  r,.,  for  the  reflection  from  the  donut  on  the  —  y-axis,  and 
for  the  reflection  from  the  donut  on  the  +  y-axis  r,;  =  fj.,.  In  Equation  (344)  and  in 
the  radius  of  curvature,  jRi,  —  and  +  signs  refer  to  reflection  from  the  donut  on  the 
—  and  -f  y-axis,  respectively.  Note  that  a  first  order  caustic  occurs  when  ^  ♦  0  and 
P2  —*  00.  This  caustic  correction  may  be  done  much  the  same  way  as  the  first  order 
diffracted  caustic  is  corrected,  which  is  discussed  in  detail  in  Section  4. 


5  Caustic  Correction  to  the  First  Order  Solution 


The  concept  of  using  equivalent  currents  for  caustic  correction  discussed  here  have 
also  been  used  in  References  [31,  32].  The  equivalent  line  currents  for  a  cylinder 
are  obtained  by  setting  the  infinite  cylinder’s  GO  reflected  field  equal  to  the  field 
radiated  by  the  infinite  line  currents.  The  equivalent  line  currents  for  a  cylinder  are 


and 


r  (345) 

/”  =  -/j  •  (»'  (0)  ■  W)  (346) 
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where 


^  a  f  -  hr 
2sin^l3' 

a  is  the  radius,  and  is  the  normal  to  the  2-D  cylinder  as  shown  in  Figure  84.  First, 
the  caustic  correction  from  half  of  the  donut  on  the  +  y-axis  will  be  considered.  The 
result  from  this  caustic  correction  wil*  be  used  to  get  the  result  from  the  whole 
donut.  Referring  to  Figure  84, 


-  ^  (d  .  A  .  . 

<‘c  —  I  -  +  a  sin  u  I  cos  <f>x  l-4-asinujsin^y  +  o  cos  9  i. 


and 


Let 


fir  =  cos  6  z  +  sin  9  y, 
$  =  sm<f)x  —  cos  ^  y, 
f  =  sin  ^  y  +  cos  9  z, 

f'  —  sin  9'  y  cos  9'  z. 


E'  =  E„x 

and  J5„  =  1.  The  far- zone  fields,  due  to  the  equivalent  line  currents,  are 

E,  =  /jjyfcsin/SV—  r  (^) 

47rr 

and 

E„^(0x  f)  jfesin/3^  fj  r  {<f,) 

Equations  (354)  and  (355)  may  be  simplified  to 

Er  =  - - sin^ 

r  27r  V  2 


and 


E„  = 


k  fa  , 


r  2ir 


gjf  ®  £  COS^  ^e>fc(«/+2o»ine).infliiin*^^ 


(348) 

(349) 

(350) 

(351) 

(352) 

(353) 

(354) 

(355) 

(356) 

(357) 
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Combining  Equations  (356)  and  (357)  and  after  some  manipulation,  it  may  be  shown 
that  the  first  order  reflected  field  with  the  caustic  correction,  is  given  as 


=  E’  (0)  , 


(358) 


where 


(359) 

(360) 


and 


=  fc  ±  asin  fir.  (f  +  f') .  (361) 

In  Equation  (360),  Jo  and  refer  to  the  Bessel  and  Hankel  functions,  respectively, 
and  in  Equation  (361),  the  +  and  the  —  signs  refer  to  the  reflection  from  donut  of 
the  +  and  —  sides  of  the  y-axis,  respectively.  The  normal  at  the  point  of  reflection 
as  defined  in  Equation  (349),  is  hr.  The  transition  function  is  Tr  {U).  The  following 
are  the  small  and  large  argument  forms  which  are  obtained  using  the  series  and  the 
asymptotic  forms,  respectively,  of  the  Bessel  and  Hankel  functions,  so  that  for 

|f^l«l, 


Tr{U) 

and  for  (f7|  >>  1, 

Tr{U) 


irU 


(362) 
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SU  128t/2 


_ fl 
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J_1 

t/2 


(363) 


6  Second  Order  Solution  and  its  Caustic  Correc¬ 
tion 


The  second  order  solution  to  the  donut  (Figure  85)  obtained  by  modification  to  the 
double  bounce  on  the  two-dimensional  cyUnder-cylinder,  discussed  in  Section  3,  is 
given  as 


=  (i’(0)  Rt)  R2 


P\ 


P2 


\  {pV  +  «)(P2  +  «) 


(364) 
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where  Ri  and  R2  are  as  defined  in  Equal  ion  (1)  at  the  first  and  second  reflection 
points,  p\'  and  p\^  are  as  defined  in  Equations  (340)  and  (341),  and 

r,  ^  (acos/3~0 

2sin  (0  + /(?)cos/9’ 

and 

1  1  2cos  (^  + /3)sin/0 

P'2  ■*  +  P2  (a  sin  /?  —  f ) 

are  obtained  using  Equation  (2). 

It  is  observed  that  as  ^  >  0,  p^  —*  00,  Equation  (364)  appears  to  become 
infinite.  This  is  expected  since  we  are  at  a  double  reflection  caustic.  However  a 
careful  evaluation  of  Equation  (364)  at  the  caustic  shows  that  the  double  reflected 
field  vanishes  at  the  caustic.  One  simple  way  to  see  this  is  to  look  at  two  double 
reflected  terms  that  are  ~  90°  apart  on  the  donut.  It  is  observed  that  along  the 
z-axis,  these  two  terms  cancel  each  other.  Thus,  the  total  resultant  field  from  the 
double  reflected  caustic,  along  the  z-axis,  vanishes.  However,  near  the  caustic,  ^  0, 

it  is  necessary  to  set  up  equivalent  currents  on  half  rims  that  are  associated  with 
the  two  reflection  points  and  integrating  the  equivalent  currents  along  half-circles, 
the  radius  of  which  are  determined  by  the  two  reflection  points  (GO),  in  order  to 
obtain  the  doubly  reflected  field.  This  interaction  could  not  be  written  in  closed 
form  because  of  the  complexity  of  the  distance  parameter,  s,  the  distance  between 
the  two  reflection  points.  (See  Figure  85.)  Performing  this  integration  numerically, 
it  was  found  that  the  doubly  reflected  term  hardly  made  any  contribution  to  the 
total  result,  at  least  for  the  dimensions  of  the  donut  considered.  (See  Figures  86  and 
87.) 

7  Results 

Figures  86  and  87  compare  the  soft  and  the  hard  backscatter  at  10  GHz  from  a  donut 
of  solid  radius,  a  =  -^  in  (0.943  X)  and  the  ring  radius,  |  (2.7  A).  Since 
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the  result  for  the  soft  polanz.ation,  Figure  86,  is  hotter  than  for  the  hard  polarisation, 
Figure  87,  it  is  expected  that  the  creeping  waves,  which  have  not  been  considered  in 
this  solution,  are  important  for  this  dimension  of  the  donut.  Since  the  result  worsen 
as  6  increases,  it  may  be  expected  that  double  bounce  mechanisms  that  are  not  in 
the  principal  plane  become  more  important. 
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SECTION  9 


UTD  Scattering  from 
Multi- Joined  Cylinders 

1  Introduction 

The  multi-joined  cylinder  model  for  the  scattering  of  the  electromagnetic  fields  is 
developed  for  the  purpose  of  studying  radar  remote  sensing  of  forested  areas.  The 
top-hat,  discussed  in  Section  5,  can  be  used  to  study  the  tree  trunk  and  ground 
interaction.  In  this  section,  the  multi-joined  cylinders  model,  see  Figure  88,  is  used 
to  investigate  the  scattering  from  tree  trunks  and  branches.  Concentric  dielectric 
layers  are  used  to  represent  the  material  characteristic  of  the  trunks  and  branches, 
and  their  shape  is  represented  by  using  more  than  one  cylinder  to  create  bends. 

Uniform  Theory  of  Diffraction  (UTD)  is  employed  to  obtain  the  scattered  fields. 
The  dihedral  study  done  in  Reference  (13]  and  outlined  in  Section  3  is  used  as  the 
building  block  for  modeling  the  cylinder-cylinder  interactions.  The  main  differences 
being  that  the  reflection  coefficients  are  modified  to  take  into  account  the  dielec¬ 
tric  layers,  see  Appendix  B  or  Reference  [33],  and  the  standard  iwo  dimensional 
UTD  diffraction  coefficients  are  modified  to  take  into  account  the  spread  factors  and 
caustic  corrections.  The  caustic  corrections  are  basically  the  same  as  discussed  in 
References  [21,  22]  and  summarized  in  Section  4.  Also,  junction  characteristics,  not 
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discussed  in  the  iiteratnre,  such  as  radius  nf  curvature  and  edge  normal,  are  required 
to  obtain  the  diffracted  fteld  from  a  junction. 


2  Theoretical  Background 

2.1  Modified  Diffraction  Coefficient 


In  order  to  find  the  scattered  fields  from  dielectric  surfaces,  the  edge  diffraction 


coefficient  for  the  perfectly-conducting  case  (Equation  (25))  discussed  in  Section  3, 
has  to  be  modified  to  take  into  account  a  junction  of  two  dielectric  plates.  It  is  shown 
in  References  (34,  35]  that  the  modified  diffraction  coefficient  for  an  electrically  thin 


dielectric  is 


— e 


2n\/2Trfc 


2n 


(367) 


where  R'J'  and  are  the  reflection  coefficient  for  the  vertical  and  horizontal 
polarization  for  the  ‘o’  and  the  ‘n’  faces,  respectively,  and  and  T*’^'  are  the 
transmission  coefficient  for  the  vertical  and  horizontal  polarization  for  the  ‘o’  and 
the  ‘n’  faces,  respectively.  Details  of  the  reflection  and  transmission  coefficient  are 
presented  in  Appendix  B  or  Reference  (33).  To  generalize  the  edge  diffraction  coef¬ 
ficient  to  include  diffraction  due  to  dielectric  plates,  the  following  assumptions  are 
made  [34] 

(a)  n  «  1,  or  n  =  2  and 

(b)  the  width  of  the  dielectric  plates  is  a  small  fraction  of  a  wavelength. 

Note  that  this  is  not  the  case  for  the  multi-cylinder  model.  However,  if  the 
cylinders  are  assumed  to  be  lossy,  the  transmission  coefficients,  T,"  *,  are  negligible. 
In  general,  for  the  backscatter  case,  the  incident  boundary  cotangent  terms  have  a 
much  smaller  contribution  as  compared  to  the  reflection  boundary  cotangent  terms. 
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Hence,  for  the  backscatter  case  the  transmission  coefficient  may  be  assumed  zero. 
Therefore,  under  the  assumption  of  backscatter  and  lossy  cylinders.  Equation  (367) 
can  be  simplified  to 


-e-n'  r 


2n\/2irk 


cot 


2n 


(368) 


2.2  The  Junction  Edge  Normal  Vector  and  Curvature 


In  order  to  find  the  diffracted  caustic  distance,  pj,  in  Equation  (22),  from  a  typical 
cylinder-cylinder  junction,  it  is  necessary  to  determine  the  edge  normal  vector,  hr, 
and  the  radius  of  curvature,  a^.  Here,  the  right  hand  x-y-z  coordinate  system  is 
defined  so  that  the  x-axis  is  along  the  normal  to  the  end  cap  of  the  main  cylinder 
(trunk)  and  the  y-axis  is  in  the  plane  formed  by  the  end  cap  normals  (also  the 
pattern  plane)  in  the  direction  of  the  junction  of  interest.  Figure  89.  It  is  shown  in 
Appendix  C  that  for  a  cylinder- cylinder  junction  angle  a  (the  dihedral  angle)  and 
equal  radii  of  cylinders,  a,  the  junction  edge  normal  is 


hr  =  cos 

and  the  radius  of  curvature  is 


o., 


(369) 


(370) 


3  Modification  to  the  Two-Dimensional  Dihe¬ 
dral 


The  only  caustic  that  needs  to  be  considered  is  that  of  single  diffractions  from  the 
end  caps.  Thus,  its  correction  is  the  same  as  that  of  the  cylindrical  part  of  the  top- 
hat,  which  has  been  discussed  in  Section  3.1.  Note  that  even  if  the  junction  angle, 
a,  is  90°  the  diffracted-reflected  and  the  diffracted-reflected-diffracted  terms  do  not 
exhibit  caustics  at  broadside. 
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Figure  89:  Cylinder-cylinder  difTracted-reflected  terms. 

3.1  Spread  Factor  and  Caustic  Correction  for  Diffracted 
Terms 

Equation  (22)  is  used  to  get  the  spread  factor  for  single  diffractions.  However,  if  the 
single  diffraction  is  from  a  junction,  nr  and  used  in  Equation  (22)  are  defined  by 
Equations  (369)  and  (370),  respectively. 

3.2  Spread  Factor  for  Diffracted-Reflected  Type  Terms 

For  the  backscatter  case,  the  spread  factor  due  to  the  diffracted- reflected  (DR)  and 
reflected-diffracted  (RD)  terms  are  the  same  since  the  two  ray  paths  are  the  same. 
For  details,  see  Section  2. 
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The  spread  factor  for  the  R  >  term  (or  RjPt )»  s****  Figure  89,  may  he  represented 


Psprcati  /actor  j  ^  ^  ^  J Fl  P2  |  » 


where  p’^,  using  Equation  (22),  is 


H  ^  ^ _ 

sin  6  +  sin  (2a  —  ’ 

and  the  principal  radii  of  curvature  of  the  reflected  wavefront  are 


p,  =  s. 


(371) 


(372) 


(373) 


1  1  2  sin  (a  —  0) 

p2  s-\-  p^  a 

or 

a(5  +  /) 

^  a  4-  (a  +  p‘^)  2  sin  (a  —  0) 
Using  Equations  (372)  to  (375)  in  (371), 


r.Hj/rrari  fnrtor 


(374) 


(375) 


(376) 


f'.prcH  fncl,^  [fl  +  (s  + /)  2  sin  (o  -  ^)  J  '  ^ 

The  spread  factor  for  the  D2R\  term  (or  Rf  D2),  see  Figure  89,  may  be  represented 


P,prrnrf  factor  [  3  (prf  +  j)  FiJ  > 

where  p’^,  using  Equation  (22),  is 

rf  ^ _ ^ _ 

sin  (a  +  5)  +  sin  (a  —  0)  ’ 

and  the  principal  radii  of  curvature  of  the  reflected  wavefront  are 

Fj  =  ■*, 


(377) 


(378) 


(379) 


a  (s  +  p'*) 
a  +  (s  +  p'')  28in  ^ 

Using  Equations  (378)  to  (380)  in  (377) 


Jhii, 

r  Mprrad  jnt  ior 


a  -t-  (.s  +  p'')2sin  9 


(380) 


(381) 
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3.3  Spread  Factor  for  Diffracted-Reflected-Diffracted  Terms 


The  spread  factor  for  the  D-iRxD^  term,  see  Figure  90,  may  be  represented  as 


ruprr.atl  factor 


pi 


P\P2  A 
P2 


[s(pf+s)  ip\-^-^){p^2  +  »y  \ 


where  pf,  using  Equation  (22),  is 

d  ^ _ a _ 

*  cosa  +  |sin(a  —  d)|’ 

and  the  principal  radii  of  curvature  of  the  reflected  wavefront  are 


and 


p\  =  s, 


1  1  2 
- rf  +  “» 


P2  «  +  Pi  ® 
and  for  the  second  diffraction,  p^,  using  Equation  (22),  is 

j  _  1  cos  a  +  |sin  (o  —  P)| 

(P2  +  a)  ^  o 


(382) 


(383) 


(384) 


(385) 


(386) 


The  spread  factor  for  the  D\R2D\  term,  see  Figure  90,  may  be  represented  as 

(387) 


Jhihih  _ 

r  nprraft  fartitr 


Pi 


P1P2 

Pi 


[3  (pf  +  3]  i.p\  +  >)  {pi  +  ^) 


where  pi,  using  Equation  (22),  is 


pi  = 


cos  Q  4-  sin  ^  ’ 

and  the  principal  radii  of  curvature  of  the  reflected  wavefront  are 


and 


Pi  = 


1  1  2 

- - ^5  +  “’ 


P2  ^  +  Pi  ® 
and  for  the  second  diffraction,  p^,  using  Equation  (22),  is 

j  1  cos  a  +  sin  0 

”  (P2  +  ^  o 


(388) 


(389) 


(390) 


(391) 
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Figure  90:  Cylinder-cylinder  DRD  terms 


4  Model  Types  and  Results 

Two  types  of  models  are  considered,  as  shown  in  Figures  91  and  92.  The  models 
are  made  of  PVC  pipes  with  absorber  placed  inside  because  it  is  inexpensive  and 
easy  to  piece  together  into  various  configurations.  They  are  chosen  to  validate  the 
dielectric  aspects  of  the  code.  The  dielectric  properties  of  the  PVC  and  the  absorber 
are  well  known.  The  configurations  are  chosen  to  resemble  the  crooks  of  trees.  It  is 
not  suggested  that  the  results  from  these  models  are  directly  comparable  to  those 
from  trees. 

To  determine  the  effects  of  terminations  while  comparing  the  Uniform  Theory  of 
Diffraction  and  the  measurement  results,  the  following  order  of  analysis  is  carried 
out.  All  the  results  are  for  the  backscatter  case.  First,  the  plastic  cylinders  and  the 
end  caps  are  covered  with  aluminum  foil;  Figures  93  and  94  show  the  vertically  and 
horizontally  polarized  fields  for  the  T-junction  at  2  GHz,  and  Figures  95  and  96  show 
the  vertically  and  horizontaUy  polarized  fields  for  the  Y-junction  at  2  GHz.  Second, 
the  plastic  cylinders  are  covered  in  aluminum  foil,  but  flat  absorbers  are  inserted 
in  the  end  caps:  Figures  97  and  98  show  the  vertically  and  horizontally  polarized 
fields  for  the  T-junction  at  2  GHz,  and  Figures  99  and  100  show  the  vertically  and 
horizontaUy  polarized  fields  for  the  Y-junction  at  2GHz.  Next,  the  cylinders  are  the 
same  as  above,  but  cone  pyramid  absorbers  are  inserted  in  the  end  taps;  Figures  101 
and  102  show  the  verticaUy  and  horizontaUy  polarized  fields  for  the  T-junction  at 
2  GHz,  and  Figures  103  and  104  show  the  verticaUy  and  horizontaUy  polarized  fields 
for  the  Y-junction  at  2  GHz. 

4.1  Specular  from  PVC  Cylinder  with  Absorber 

Figure  105  shows  the  real  and  the  imaginary  parts  of  the  dielectric  constants  for  the 
absorber  used  to  fiU  the  PVC  pipe.  The  dielectric  constant  of  the  PVC  pipe  is  taken 
to  be  that  of  Polystyrene  given  as  e  =  2.55  —  j%.5E  —  04  [36]  and  is  assumed  fixed 
over  2-18  GHz. 


150 


Figure  91:  The  T-junction  model. 


Figure  92:  The  Y- junction  model. 
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Figure  93:  Backscatter  from  T-junction  covered  in  aluminum  foil  at  2  GHZ  for 
vertically  polarized  field. 


Figure  94:  Backscatter  from  T-junction  covered  in  aluminum  foil  at  2  GHZ  for 
horizontally  polarized  field. 
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Figure  95:  Backscatter  from  Y-junction  covered  in  aluminum  foil  at  2  GHZ  for 
vertically  polarized  field. 


Figure  96:  Backscatter  from  Y-junction  covered  in  aluminum  foil  at  2  GHZ  for 
horizontally  polarized  field. 
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Figure  97:  Backscatter  from  T-junction  covered  in  aluminum  foil,  but  flat  absorbers 
inserted  in  end  caps,  at  2  GHZ  for  vertically  polarized  field. 


Figure  98:  Backscatter  from  T-junction  covered  in  aluminum  foil,  but  flat  absorbers 
inserted  in  end  caps,  at  2  GHZ  for  horizontally  polarized  field. 


154 


8 


vro  soiwnoM 
ICASMOOn 


Figure  99:  Backseat  ter  from  Y-junction  covered  in  aluminum  foil,  but  flat  absorbers 
inserted  in  end  caps,  at  2  GHZ  for  vertically  polarized  field. 


Figure  100:  Backscatter  from  Y-junction  covered  in  aluminum  foil,  but  flat  absorbers 
inserted  in  end  caps,  at  2  GHZ  for  horizontally  polarized  field. 
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Figure  101:  Backscatter  from  T-junction  covered  in  aluminum  foil,  but  cone  pyramid 
absorbers  inserted  in  end  caps,  at  2  GHZ  for  vertically  polarized  field. 


Figure  102:  Backscatter  from  T-junction  covered  in  aluminum  foil,  but  cone  pyramid 
absorbers  inserted  in  end  caps,  at  2  GHZ  for  horizontally  polarized  field. 
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Figure  104:  Backscatter  from  Y-junction  covered  in  aluminum  foil,  but  cone  pyramid 
absorbers  inserted  in  end  caps,  at  2  GHZ  for  horizontally  polarized  field. 
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In  order  to  match  the  theory  with  the  measurements  in  the  specular  region  broad¬ 
side  to  the  PVC  pipe  with  absorber  flat  at  the  end  caps  it  is  necessary  to  set  an  air 
gap  between  the  PVC  and  absorber  in  both  the  UTD  and  the  eigenfunction  solution. 
This  may  be  explained  by  the  fact  that  the  absorber  is  not  tightly  fitted  into  the 
PVC  pipe.  Figures  106  and  107  show  the  horizontal  and  the  vertical  polarizations 
from  a  24  in  long  cylindrical  pipe  of  radius  1.5  tn.  An  air  gap  of  0.18  m  is  chosen 
for  the  theory  to  match  the  measured  results.  It  is  observed  that  the  PVC  pipe 
was  slightly  elliptical  which  can  explain  why  the  ripple  in  the  eigenfunction  solu¬ 
tion  due  to  the  creeping  wave  does  not  match  with  the  creeping  wave  ripple  in  the 
measurement. 


4.2  Results  of  PVC  pipe  with  absorber 


In  the  modified  diffraction  coefficient,  Equation  (368),  one  face  of  the  wedge  is 
unaware  of  the  characteristics  of  the  other  face.  This  suggests  that  simplifying 
Equation  (368)  to  retain  information  of  only  one  face  at  a  particular  direction  should 
provide  even  better  results.  The  diffraction  coefficient  is  simplified  as 
if  cot  I  >  Ik-  cot  |, 


D,.h  = 


—  1 


2n\/2irk  " 


cot 


/'tt  -  (# -^  $')\ 

V  2n  )’ 


D,,h  = 


2ny/2irk  \  2n  } 


(392) 


2n\/2'nk 

It  is  found  that  patterns  and  signatures  of  cylinders,  T  and  Y  junctions,  ob¬ 
tained  by  using  the  simplified  diffraction  coefficient.  Equation  (392),  compare  better 
with  measurements  than  those  obtained  by  using  the  modified  diffraction  coefficient, 
Equation  (368).  For  all  the  results  shown  in  this  subsection,  the  simplified  diffraction 
coefficient  of  Equation  (392)  is  used  and  an  air  gap  of  0.18  in  is  assumed  between 
the  PVC  pipe  and  the  absorber.  Also,  all  the  results  are  for  the  backscatter  case  at 
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Figure  106:  Specular  from  24  in  long  cylinder  of  radius  1.5  in  for  horizontally 
polarized  field. 


Figure  107:  Specular  from  24  in  long  cylinder  of  radius  1.5  tn  for  vertically  polarized 
field. 
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2  CTTz.  For  a  21  in  long  r3’lindriral  pipe  of  radius  1.5  in,  Figures  108  and  100  show 
the  vertical  and  the  horizontal  polarizations,  respectively.  Figures  1 10  and  111  show 
the  vertical  and  the  horizontal  polarizations,  respectively,  for  a  cylindrical  pipe  of 
the  same  length  but  of  radius  2  in.  For  the  T-junction  model,  shown  in  Figure  91, 
Figures  112  and  113  show  the  vertical  and  the  horizontsd  polarizations,  respectively. 
Figures  114  and  115  show  the  vertical  and  the  horizontal  polarizations,  respectively, 
for  the  Y-junction  model,  shown  in  Figure  92. 

5  More  Realistic  Modeling 

The  PVC  models  studied  above  have  been  chosen  simply  to  validate  the  theory  and 
the  computer  code.  More  work  needs  to  be  done  to  come  up  with  a  reasonable  model 
for  a  tree.  Certain  conclusions  from  this  study,  however,  may  be  pointed  out.  Since 
the  PVC  pipe  is  not  lossy,  it  is  expected  that  better  results  can  be  obtained  for 
actual  trunks  and  branches.  This  is  because  the  surface  waves  will  dampen  out  due 
the  wood  being  lossy.  Also,  in  general,  it  is  expected,  that  the  truncations  of  the 
cylinders  which  are  considered  in  this  study  only  approximately  model  the  actual 
shapes  of  the  trunks  and  branches.  One  needs  to  come  up  with  an  equivalent  length 
in  order  to  model  a  trunk  or  a  branch. 
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Figure  109:  Backscatter  from  24  in  cylinder  of  radius  1.5  in  at  2  GHZ  for  horizontally 
polarized  field. 


Figure  110:  Backscatter  from  24  in  cylinder  of  radius  2.0  in  at  2  GHZ  for  vertically 
polarized  field. 
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Figure  114:  Backscatter  from  Y-junction  at  2  GHZ  for  vertically  polarized  field. 


Figure  116:  Backscatter  from  Y-junction  at  2  GHZ  for  horizontally  polarized  field. 
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SECTION  10 


Conclusion 


In  general,  the  high  frequency  EM  scattering  from  complicated  physical  structures 
such  as  aircraft,  tanks,  trucks  etc.  can  be  analyzed  using  ray  techniques  by  building 
up  such  complex  structures  from  simpler  shapes  such  as  cylinders,  cones,  dihedrals 
etc..  Therefore  a  complex  problem  may  be  analysed  piece  by  piece  and  then  applying 
superposition.  In  this  work,  classical  Uniform  Theory  of  Diffraction  (UTD)  with 
modifications  for  the  caustics  and  curvatures  is  used  to  gel  the  scattered  fields. 
Shapes  where  the  second  and  higher  order  interactions  are  important  over  a  wide 
range  of  angle  are  considered. 

It  is  shown  that  for  the  rectangular  dihedral  in  the  principal  plane,  the  order  of 
interactions  required  to  get  the  pattern  in  the  dihedral  region,  for  a  dihedral  of  angle 
a,  is  found  to  be  the  integral  value  of  for  the  backscatter  case.  Since  the 

computer  code  developed  for  this  analysis  contains  all  interactions  up  to  third  order 
except  the  triple  diffractions,  the  solution  discussed  here  is  valid  for  backscatter  for 
dihedrals  of  a  —  60°  or  greater.  This  solution  is  expected  to  degrade  gracefully 
in  the  dihedral  region  for  a  less  than  60°.  In  general,  as  the  angle  of  the  dihedral 
diverts  from  a  which  makes  (^~)  integer,  the  main  beam  amplitude  falls  fairly 
rapidly.  The  UTD  formulation  provides  excellent  results  compared  with  the  method 
of  moments.  For  the  case  of  maximum  scattering,  from  a  GO  point  of  view,  and 
for  a  general  angled  dihedral  and  bistatic  angle,  the  dominant  UTD  terms  in  the 
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dihedral  region  are  reduced  t.o  simple  forms.  Comparisons  of  these  dominant  tlTD 
results  with  the  total  UTD  solutions  show  good  agreement. 

The  two-dimensional  result  of  the  rectangular  dihedral  is  modified  by  caustic 
corrections  and  spread  factors,  to  get  three-dimensional  results  for  the  top-hat,  bi¬ 
cone,  and  parallel  plate  and  circular  cavities.  It  is  found  that  the  response  of  cavities 
is  not  as  sharp  as  compared  to  that  of  the  90°  dihedral.  The  results  compare  well 
with  the  those  obtained  from  body  of  revolution  moment  method  solution. 

Scattering  from  the  donut  is  considered  in  Section  8.  The  first  and  second  order 
GO  field  caustics  are  corrected  by  using  equivalent  currents.  Again,  the  results 
compare  well  with  the  body  of  revolution  moment  method  solution. 

Scattering  from  cylindrical  pipes  that  form  a  T  and  Y-junctions  of  perfect  electric 
conductor  and  dielectric  materials  are  considered.  It  is  shown  that  by  using  the 
modified  diffraction  coefficients,  the  UTD  results  compare  well  with  measurements 
and  an  eigenfunction  solution.  Some  more  work  is  required  to  extend  these  results 
for  scattering  from  tree  trunks  and  branches. 

It  is  observed  that  a  few  underlying  principles  are  all  that  is  required  to  get 
scattered  fields  from  many  different  shapes.  It  is  expected  that  even  better  agree¬ 
ment  and/or  results  for  small-sized  models  could  be  achieved  by  including  more 
interactions. 
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Appendix  A 

Balancing  of  Fields  for  Joined 
Wedges 


In  this  appendix,  a  detailed  mathematical  analysis  is  carried  out  in  order  to  obtain 
various  diffracted  terms  which  are  then  used  in  Section  4  to  show  continuity  of  the 
fields.  Reference  is  made  to  Figure  27.  For  diffraction  from  Qj  the  angles  are 


</?'  =  TT  +  et  ond  —  ^2- 


(393) 


Thus  diffraction  from  Q\  is  given  as 

r  2n 


,  (!^) ,  ^}] 


(394) 


Since  the  diffracted  field  is  taken  close  to  the  incident  shadow  boundary,  jPc,,n(t?i) 
may  be  split  into  a  term,  FJ„,(Qi),  which  is  singular  at  the  shadow  boundary,  and 
non  -  singular  term  that  is  denoted  by  F”,„(Qi). 


1 


1 


Lcj  -  ei  62  +  6] 

Noting  that  P,,m[Q\)  is  present  only  when  ej  >  0 


(395) 


KJQ<)  = 


e-Wi 

V2nk 


e;  -  ,  -j 


e  4  flsrn(f2)] 


(396) 
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and 


X  r  /27r  +  e, -e2^ 


27r  +  £)  —  €2' 


T  cot 


2n\/2Tik 
(2r  4-  Cl  -f 


2n 


e  2^^  +J>5n(f2)]. 


The  angle  parameters  for  Q2  are 

4>'  —  — t)  and  <l>  =  IT  —  62 

and  the  diffraction  from  Q2  is 

-g-W  r  /2ir-C2  +  ci 


(397) 


(398) 


n2  n  rr\  \  \  l  ~  ^2  +  f|  \  , 

A,.„  =  />.  .,.(<?.)  =  [col  ( - - j  +  3 


2n 


{col( 


2n 

22r  —  €2  —  Cl 


-6|  +  £2 


)+-^y 

/  +  1 1  J  . 


Again  separating  (<?2)  in  to  singular  and  the  non-singular  parts  and  noting  that 
is  only  present  when  ei  <  0. 

ttttJ  '‘'*1  '>  -  <«•» 


,-jTr/l 


d"  \  f  ^2jr  — e2  +  ei\ 

r'  \  2n  -} 


(401) 


|i )l 

(402) 

The  angle  parameter  for  Qj  are 

1 

(!>'  =  7r/2  4-  £i 

and  <!>  =  37r/2  -1-  t> 

(403) 

For  the  single  diffracted  field  from  Q,i 

.  Note  ONLY  the  term  that 

is  singular  for 

diffraction  form  Q:\  is  shown  below.  Thi.s  term  is  associated  with  ir  -  (^  - 


p:jq<)  = 


s/Tnk 


f,  -  C2 


’kH 


(404) 
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From  Figure  27  and  Equation  (391),  we  note  that  for  the  diffrartion  from  edge  Qi , 
the  term  associated  with  tt  —  —  blows  up.  From  double  diffraction  formulation, 

n  —  is  obtained  by  setting  p  =  1  in  Equation  (42).  Again  from  Figure  27  and 

Equation  (399),  it  is  noted  that  for  diffraction  from  edge  Q2  the  term  associated  with 
IT  +  (^  —  <l)')  blows  up.  Note  that  in  the  double  diffraction  formulation,  the  angles 
<i>  and  <!>'  on  the  second  edge  {Q2)  are  defined  in  the  reverse  to  those  defined  for 
single  diffraction.  With  this  in  mind,  for  double  diffraction,  ir  ~  {(f>  —  <f>')  is  obtmned 
by  setting  g  =  1  in  Equation  (43).  Thus,  we  expect  the  double  diffraction  term 
obtained  by  setting  p  =  1  and  g  =  1  in  Equation  (44)  to  compensate  for  the 
discontinuity  and  singularities  of  the  single  diffracted  fields. 

For  the  double  diffraction  formulation 


and 


Thus 


<i>\  =  (ni  -  l)ir  -  e, 

(405) 

^'2  -  (”2  “  1)’^  +  «2- 

(406) 

^p=i  =  TT  -  (niTT  -  (n,  -  l)ir  +  Ci]  =  -€| 

(407) 

=  TT  -  [rijTT  -  (n2  -  l)7r  -  62)  = 

(408) 

ap=i  =  2cos^{— £1/2  —  7r/2]  =  25m^(ei/2) 

(409) 

aq=i  =  2co5‘^[c2/2  —  7r/2)  =  2sm^(£2/2) 

(410) 

Using  the  above  equations  in  Equation  (46) 

F.(V.,$,=  ,)  =  (— )  (— ) 

47rjnin2  \  Ci  /  \  €2  J 

2sin"(c2/2)F[ikd2sin"(c,/2)]  -  2sin''(e,/2)Fllfed2sin"(e2/2)] 

2sm^(£,/2)  -2sin^(e2/2) 

\kVd)- 


For  small  x 


F[x]  -  [y/^  -  2a  e-"/') 


(411) 


(412) 
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F[2fc<ism^(ei /2)]  = 


^2Tfcd(e,/2)2  -  4ifc(/(e,/2)^e^"/'] 

1  e,  I 


~j(^l  l+2W(f/l) 
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(413) 

(414) 
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neglecting  terms  of  order  Cj  and  C2 
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From  Equation  (45) 

/^.(</ip=,,4i,=  .)  =  (^1^) 
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t\l2~4l2  [ky/dj 


(418) 


Using  Equation  (414) 
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neglecting  terms  of  order  £1  and  £2 
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Appendix  B 

Reflection  and  Transmission  of 
Plane  Waves  in  Planar,  Stratifled 
(Multi-Layered)  Medium 


Details  of  reflection  and  transmission  of  the  plane  wave  in  multi-layered  medium  are 
presented  in  (33).  Here  a  summary  of  (33)  is  presented.  Suppose  a  plane  wave  in 
free  space  has  oblique  incidence  on  a  plane  multilayer  consisting  of  N  homogeneous 
isotropic  slabs  as  shown  in  Figure  116.  Let  d„,  and  e  represent  the  thickness, 
permeability,  and  permittivity  of  slab  n.  The  slabs  are  considered  to  have  infinite 
width  and  hight  and  parallel  surfaces,  with  unbounded  free  space  on  both  sides  of 
the  multilayer.  The  incident  plane  wave  impi'iging  on  the  left-hand  surface  of  the 
multilayer  is  given,  in  the  TE  case  (i.e.  perpendicular  polarization)  by 

E;  =  (421) 

where  6  is  the  angle  of  incidence,  fco  =  2irf\,  and  A  is  the  free  space  wavelength. 
The  reflected  wave  is  given  by 

E;  =  (422) 

where  R  is  the  reflection  coefficient  of  the  multilayer.  The  transmitted  plane  wave 
on  the  right-hand  side  of  the  multilayer  is  represented  by 

(423) 
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Figure  116:  A  plane  multilayer,  illustrating  the  outgoing  and  the  reflected  waves  in 
each  layer. 

where  T  is  the  transmission  coefficient  of  the  multilayer.  The  field  in  each  layer  can 
be  regarded  as  an  infinite  series  of  the  plane  waves  bouncing  back  and  forth,  but  it 
is  more  convenient  (and  equally  valid)  to  consider  it  to  be  the  sum  of  only  two  plane 
waves,  one  travelling  outward  and  one  reflected.  In  layer  n,  for  example,  the  field  is 
represented  by 

E„  =  e-’""*'"”".  (424) 

Similarly,  in  layer  n+1  the  electric  field  intensity  is  given  by 

E,.^,  =  (425) 

The  boundary  between  layers  n  and  n  +  1  is  located  at 

z„  =  dt  -i-  d:\  +  -) — h  d„.  (^26) 


By  enforcing  the  boundary  conditions  on  E^  and  Hy  at  z  =  2„,  it  can  show  that 
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Q„ ' 

’  A„  ' 
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where 
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“  2 

( 1  1 

1  ~'>n+ 1  )*n 

(428) 
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(  M»7"+i  / 

1  7n4  I  )■*" 

1  ^  ’ 

(429) 
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1 

”  2 

/  J  /^n+l7n 

\  /^n7n+1  , 

1  *('rn+7n+i)^n 

(430) 

s„ 

j  1  A*n+l7t)\ 
/^n7ri  +  l/ 

g"  ('Vti  “'>n+  1  )^n 

(431) 

The  propagation  constant  7,,  for  layer  n  will  be  complex  if  the  medium  is  dissi¬ 
pative.  Both  the  real  and  the  imaginary  parts  of  7,,  will  be  positive.  If  layer  n  is 
a  lossless  medium,  7,,  will  be  purely  imaginary.  The  wave  equation  is  employed  to 
obtain 

7fi  =  (432) 

The  reflection  and  transmission  coefficient  of  the  multilayer  can  be  calculated  in 
a  systematic  manner  by  setting 

A. =  1  (433) 

and 

Bo  =  0  (434) 

and  then  by  using  the  recursion  equations  Equation  (427)  to  calculate  >1),  Bi,  A2, 
Bj,  ...  and  Ba  +  i  in  that  order. 

From  Eqs.  (421)  through  (424), 


and 


Eu  —  j4a+i> 

(435) 

(436) 

R=-j^- 
■^A  +1 

(437) 
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In  the  TE  case  the  constants  A„  and  B„  represent  the  electric  field  intensities  of  the 
outgoing  and  reflected  wave  in  each  layer.  In  the  TM  case  (parallel  polarization) 
the  solution  proceeds  in  the  same  manner.  The  equations  given  above  apply  in  both 
cases  but  the  A„  and  B„  represent  the  magnetic  field  intensities  in  the  TM  case  and 
/i„+i  and  p„  must  be  replaced  with  e„^.|  and  €„  in  Eqs.  (428)  through  (431). 

If  a  perfectly  conducting  sheet  is  placed  on  the  right-hand  surface  of  the  multi¬ 
layer  (i.e.,  on  the  x-y  plane),  the  solution  is  again  given  by  the  equations  above  with 
the  exception  that  the  transmission  coefficient  T  is  not  calculated  in  this  case,  and 
Eqs.  (433)  and  (434)  are  replaced  with 

A)  =  1  and  Bt  =  —1  tn  the  TE  case  (438) 

and 

A)  =  1  and  B|  =  1  in  the  TM  case.  (439) 

Equations  (438)  and  (439)  are  obtained  by  forcing  the  tangential  electric  field  inten¬ 
sity  to  vanish  at  the  perfectly  conducting  plane. 

In  the  above  equations,  the  reflection  coefficient  R  is  defined  as  the  ratio  of  the 
reflected  wave  amplitude  to  the  incident  wave  amplitude  at  the  coordinate  origin; 
that  is, 

R  =  M  the  TE  case,  (440) 

Bj,  (0,0,0)’  ^  ' 

and 

R  =  /<”■  the  TM  case.  (441) 

(0,0,0)  ^  ^  ’ 
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Appendix  C 

Cylinder- Cylinder  Junction  Edge 
Normal  and  Curvature 


Consider  a  “natural”  cylinder-cylinder  junction  shown  in  Figure  117.  The  cylinders 
have  radii  a/  and  02  and  make  an  angle  a  at  the  junction  of  interest.  The  position 
vector  of  cylinder  1  in  cylindrical  coordinates  is 


fi  =  a\  (cos  •+  sin  y, )  4-  2i ii 


(442) 


and  of  cylinder  2 

f2  =  02  (cos  <^2*2  +  sin  (^2^2)  +  2222>  (443) 

where  <f>\  is  measured  from  ii-axis  and  (f>2  is  measured  from  the  a;2-axis  in  the  X\  —  yi 
and  X'i  — 1/2  planes,  respectively.  The  two  right  handed  coordinate  systems  {xf,yf,Zi 
and  x-iiyziZ-z)  are  related  so  that 


X2  —  —  cos  axi  —  sin  azi , 


(444) 


Z2  =  sin  ail  —  cos  aij 


(445) 


and 


y?  =  Vi- 


(446) 


Using  Equations  (444)-(446)  in  Equation  (443) 


r2  =  (22  sin  a  —  a-i  cos  <^2  cos  a)  i|  -|-02sin  (l>  ,y\  -  (02  cos  <f>2  sin  a  -t-  22  cos  a)  2| ,  (447) 
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and  at  the  junction 

f,  =  fj. 

(448) 

Now,  the  junction  may  be  defined  by: 

a\  cos  <^|  —  z-i  sin  a  —  a-z  cos  4>z  cos  a, 

(449) 

Uj  sin  <j^\  —  0,2  sin  (j>2 

(450) 

and 

2i  =  —  (a2  cos  <f)2  sin  a  4  zi  cos  a) . 

(451) 

r 

From  Equations  (449)  and  (450)  one  gets 

0  =  zl  sin^  a  —  22^2  cos  <(>2  sin  (2q)  4  {cos^  <f>2  cos'^  a  4  sin^  ^2)  —  Op 

(452) 

Solving  the  above  Equation  for  22, 

22  =  02  COS  <j^2  cot  a - r —  \/af  -  oj 

SI  1  a  ^ 

(453) 

181 


If  the  radii  of  the  cylinders  are  the  same,  at  —  a2  =  a,  then  Equation  (453)  may  be 
simplified  to 

acos<^-2 

22  =  — : - Icosa  —  II  (454) 

sin  a 


-  /  /  *  •  /  -  \  ®  *^OS  <p2  ,  , ,  . 

r2  =  a  (cos  <P2®2  +  sin  <P2y2)  4 - : - [cosa  —  IJ  Z2- 

sin  a 


(455) 


The  first  and  the  second  derivatives, 


Ej  =  — ^  =  a  (—  sin  4>2X2  +  cos  ^2^2)  —  ^  q  -  1)  ij  (456) 

d(f>2  sin  a 


f "  = 


/  *  •  ,  *  \  ®  ^OS  <f)2  , 

— —  =  -a  (cos  <P2X2  +  sm  (pzili) - : - (cos  q  —  1)  221  (457) 

d<p2  sin  Q 


are  needed  in  order  to  determine  the  curvature  and  the  edge  normal  vector. 
At  the  point  of  diffraction,  the  junction  radius  of  curvature  is 


I  If]''  fa\ 

'’I**’-' =  =  “">'12  j  ’ 


and  the  edge  normal  vector,  n,.,  is 


K  =  U,=0  =  cos  (I)  X2  -  sin  22 


(458) 


(459) 


n,  =  cos 


(460) 
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